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We propose a criterion for the validity of semiclassical gravity (SCG) which is based on the stability 
of the solutions of SCG with respect to quantum metric fluctuations. We pay special attention to the 
two-point quantum correlation functions for the metric perturbations, which contain both intrinsic 
and induced fluctuations. These fluctuations can be described by the Einstein-Langevin equation 
obtained in the framework of stochastic gravity. Specifically, the Einstein-Langevin equation yields 
stochastic correlation functions for the metric perturbations which agree, to leading order in the 
large A'' limit, with the quantum correlation functions of the theory of gravity interacting with N 
matter fields. The homogeneous solutions of the Einstein-Langevin equation are equivalent to the 
solutions of the perturbed semiclassical equation, which describe the evolution of the expectation 
value of the quantum metric perturbations. The information on the intrinsic fiuctuations, which 
are connected to the initial fluctuations of the metric perturbations, can also be retrieved entirely 
from the homogeneous solutions. However, the induced metric fluctuations proportional to the noise 
kernel can only be obtained from the Einstein-Langevin equation (the inhomogeneous term). These 
equations exhibit runaway solutions with exponential instabilities. A detailed discussion about 
different methods to deal with these instabilities is given. We illustrate our criterion by showing 
explicitly that flat space is stable and a description based on SCG is a valid approximation in that 
case. 

I. INTRODUCTION 

In this paper we discuss the conditions underlying the validity of semiclassical gravity (SCG) emphasizing the role 
of metric fluctuations induced by the quantum matter sources. SCG is based on the self-consistent solutions of the 
semiclassical Einstein equation for a classical spacetime driven by the expectation value of the stress tensor operator 
of quantum matter fields. We propose a criterion based on stochastic semiclassical gravity 0, 0] and compare it 
with the recently proposed criterion by Anderson et al. based on linear response theory. To do this we need to 
reexamine all the relevant factors old and new contributing to this issue, such as the reduction of higher derivative 
equations, intrinsic and induced fluctuations, and the relation between stochastic and quantum correlations. It also 
necessitates some clarification of the relation between our approach based on stochastic dynamics and the linear 
response approach, and differences with the approach pursued by Ford et al. |lt|S|, 6, 7, 8, ^ 0, [13 based on the 
normal-ordering and integration-by-parts procedures on the stress-energy bitensor. The connection clarified and the 
bridges built in this process are beneficial to further development of the "bottom-up" approaches to quantum gravity 
starting from SCG [ll 111 El El . 

a. Metric fluctuations SCG accounts for the averaged back reaction of quantum matter fields and can be regarded 
as a mean field approximation that describes the dynamics of the mean spacetime geometry. However, it does not 
account for the effects of the fluctuations of spacetime geometry, which can also be very important. Consider, for 
instance, the metric fluctuations induced by the vacuum fluctuations of the inflaton field in inflationary cosmological 
models. Those fluctuations play a crucial role in the generation of the primordial inhomogeneities which gave rise 
to the large scale structure of the present universe as well as the observed anisotropies of the cosmic microwave 
background of radiation. 

This paper focuses on the effects of the quantum fluctuations of the metric. We will restrict our treatment to small 
metric perturbations around a given background geometry. (Of course, a full treatment of those fluctuations would 
require a complete theory of quantum gravity). We will linearize and quantize those metric perturbations including 
their interaction with the quantum matter flelds. This can be described more precisely in terms of A^ identical matter 
fields. Our approach corresponds then to computing the quantum correlation functions for the metric perturbations 
to leading order in a 1/A^ expansion. 

In fact, one can show that the leading order contribution to the quantum correlation functions in a large A^ expansion 
is equivalent to the stochastic correlation functions obtained in the context of stochastic semiclassical gravity. Whereas 
SCG is based on the semiclassical Einstein equation with sources given by the expectation value of the stress tensor 
operator of the quantum matter fields, stochastic semiclassical gravity is based on the Einstein-Langevin equation, 
which has in addition sources due to the noise kernel. The noise kernel is the symmetrized connected part of the 
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two-point quantum correlation function of the stress tensor operator with respect to the state of the matter fields, 
and describes their stress-energy fluctuations. 

Making use of the equivalence between quantum and stochastic correlation functions in stochastic semiclassical 
gravity, one is naturally led to separate the symmetrized quantum correlation function for the metric perturbations 
(to leading order in into two separate contributions: the intrinsic and the induced fluctuations. The former are 

connected to the dispersion of the initial state of the metric perturbations, whereas the latter are proportional to the 
noise kernel and are induced by the quantum fluctuations of the matter fields' stress tensor operator. 

b. Validity of semiclassical gravity Different aspects concerning the validity of the description provided by SCG 
in the case of free quantum matter fields in the Minkowski vacuum state propagating on Minkowski spacetime have 
been studied by a number of authors. Most of them considered the stability of such a solution of SCG with respect to 
small perturbations of the metric. Horowitz was the first one to analyze the equations describing those perturbations, 
which involved higher order derivatives (up to fourth order), and found unstable solutions that grow exponentially 
with characteristic timescales comparable to the Planck time This was later reanalyzed by Jordan with 
similar conclusions jl^ . However, those unstable solutions were regarded as an unphysical artifact by Simon, who 
argued that they lie beyond the expected domain of validity of the theory and emphasized that only those solutions 
which resulted from truncating perturbative expansions in terms of the square of the Planck length are acceptable 
[T9I |20| . Further discussion was provided by Flanagan and Wald , who advocated the use of an order reduction 
prescription first introduced by Parker and Simon 22] but insisted that even nonperturbative solutions of the resulting 
second order equation should be regarded as acceptable. Following these approaches Minkowski is shown to be a stable 
solution of SCG with respect to small metric perturbations. 

Anderson, Molina-Pan's and Mottola have recently taken up the issue of the validity of SCG again. Their starting 
point is the fact that the semiclassical Einstein equation will fail to provide a valid description of the dynamics of the 
mean spacetime geometry whenever the higher order radiative corrections to the effective action, involving loops of 
gravitons or internal graviton propagators, become important (see Refs. 113, 13 113 for some attempts to include 
those effects). Next, they argue qualitatively that such higher order radiative corrections cannot be neglected if the 
metric fluctuations grow without bound. Finally, they propose a criterion (a necessary condition) to characterize 
the growth of the metric fluctuations, and hence the validity of SCG, based on the stability of the solutions of the 
linearized semiclassical equation. 

c. Our criterion In this paper we addresses the issue of the stability of semiclassical solutions with respect to small 
quantum corrections. When the metric perturbations are quantized, the semiclassical equation can be interpreted 
as the equation governing the evolution of the expectation value of the operator for the metric perturbations. We 
introduce a stability criterion based on whether the metric fluctuations grow without bound or not by considering the 
behavior of the quantum correlation functions of the metric perturbations. Furthermore, we emphasize that one should 
consider not only the intrinsic fluctuations, but also the induced ones. In fact, the induced fluctuations play a crucial 
role when considering the stability of simple open quantum systems for several reasons. First, those systems usually 
exhibit a characteristic relaxation time so that for much larger times the contribution from the intrinsic fluctuations 
becomes negligible. Second, after that transient period the stability around an equilibrium configuration is the result 
of a balance between the energy dissipated by the system and the fluctuations induced by the environment, which is 
encoded in the so-called fluctuation-dissipation relation connecting the dissipation and the noise kernels. 

It is true that the effect of intrinsic fluctuations can be deduced from an analysis of the solutions of the perturbed 
semiclassical Einstein equation, but in general one cannot retrieve the effect of the induced fluctuations from it. This 
effect can be properly accounted for in the stochastic semiclassical gravity framework. Both intrinsic and induced 
fluctuations are innate in the Einstein-Langevin equation. 

d. Ford's program Ford P| was amongst the first to have noted the importance of quantum fluctuations in these 
issues. An earlier criterion put forth by Kuo and Ford ^ used the variance of the fluctuations of the stress tensor 
operator compared to the mean value as a measure of the validity of SCG. As pointed out by Hu and Phillips [2^12^ 
(see reply by Ford and Wu 0) such a criterion should be refined by considering the back reaction of those fluctuations 
on the metric. Ford and collaborators also considered both intrinsic ("active") [oL ITol ITlj and induced ("passive") 
^ ^5 X fluctuations, but they did not treat them in a unifled way and did not discuss their precise relation to 
the quantum correlation function for the metric perturbations. Furthermore, they did not include the full averaged 
back reaction of the matter flelds self-consistently, and the contribution from the vacuum fluctuations in Minkowski 
space were discarded. As these issues have been discussed before by both groups of Ford and Hu, we will only make 
a few remarks at the end of this paper. 

Here, our attention will be focused on comparing the criteria based on the linear response approach proposed by 
Anderson et al. and our stochastic gravity approach. Since the differences in the two ways to address the issue of the 
validity of SCG is rooted in the difference between linear response theory and stochastic dynamics as applied to SCG, 
we hope that this work can also serve the purpose of offering a comparison between these two important approaches 
exploring the validity of the mean field approximation. In the examples provided, we will specialize the matter fields 
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to the case of free scalar fields, but generalization to vectorial or fermionic fields should not pose major difficulty. 

e. Terminology and organization To avoid unnecessary ambiguities or confusion in interpretation, it is useful to 
clarify the use of some terminology here. 

First, a comment on the difference between the stochastic gravity program in general and its present implementation 
status in particular. Stochastic semiclassical gravity can be understood as the Gaussian approximation to stochastic 
gravity. Although technically the actual implementations of stochastic gravity so far (to which our present discus- 
sion applies) have been restricted to linear metric perturbations around the background geometry and a Gaussian 
stochastic source, the theoretical construct of the stochastic gravity program has a much broader meaning beyond 
these limitations. It refers to the range of theories based on the second and higher order correlation functions. Noise 
can be defined in fully nonlinear theories {e.g. correlation noise p9l | in the Schwinger-Dyson equation hierarchy) to 
some degree |80j . but one should not expect the simple Langevin form with Gaussian and additive noise to prevail. 
Thus, stochastic gravity in this broad sense entails the whole hierarchy of correlation functions, which would imply 
going beyond order 1/N in the generating functional. It could in principle provide the means (similar to the BBGKY 
hierarchy in kinetic theory) to access the full theory of quantum gravity [T3. ITEI | . It is in this sense that we say 
stochastic gravity is the intermediate theory between SCG (a mean field theory based on the expectation value of 
the energy momentum tensor of quantum fields) and quantum gravity (understood as the full hierarchy of correlation 
functions retaining complete quantum coherence). 

Second, the precise meaning in our use of the terms perturbations and fluctuations. By perturbations of the 
metric we mean deviations of the perturbed metric from a background metric. Perturbations are purely classical and 
deterministic in general relativity and SCG. In stochastic gravity, they are classical but stochastic (with a vanishing 
statistical expectation value) so that the background configuration can be regarded as the expectation value of a 
stochastic metric (a complete gauge fixing is required to meaningfully talk about the expectation value of a metric). 
In linear quantum gravity, perturbations are quantum operators. For a state with a vanishing expectation value, the 
background metric can then be regarded as the expectation value of the metric operator (again a complete gauge 
fixing is required) times the identity operator. On the other hand, the term fluctuations is employed only to refer 
to the statistical fluctuations of the metric perturbations when they correspond to a stochastic process, or to the 
quantum fluctuations of the metric perturbations when they are treated as a quantum operator. 

Third, by leading order in the large N limit we mean the lowest order in 1/A^ with a nonvanishing contribution. 
Hence, as we will see, the leading order for the source of the semiclassical Einstein equation, which is proportional 
to the expectation value of the stress tensor operator, is whereas the leading order for the quantum two-point 

correlation functions is 1/iV. 

The paper is organized as follows. In Sec. |n]we briefly review the fundamental aspects of SCG and how one can 
study linearized perturbations around a background solution of SCG. This is generalized to incorporate the metric 
fluctuations in Sec. lIIII where the key elements of stochastic semiclassical gravity are introduced and the equivalence 
between stochastic and quantum correlation functions is explained. In Sec. IIVI we propose a generalized stability 
criterion that includes the metric fluctuations, which is then applied to the specific case of a Minkowski background. 
We conclude by summarizing and discussing the main results in Sec. 

A number of additional details and technical points are left for the Appendices. In Appendix 1X1 we illustrate the 
basic aspects of intrinsic and induced fluctuations using a simple quantum Brownian motion model. In Appendix IbI 
we provide the expressions for the dissipation and noise kernels in a Minkowski spacetime and the vacuum state. Some 
of the main steps to show the equivalence between stochastic and correlation functions using a large N expansion 
are summarized in Appendix [CI The physical interpretation of the singular coincidence limit for the noise kernel and 
possible ways to deal with it are explained in Appendix^ Finally, in Appendix^we discuss the existence of runaway 
solutions in SCG and stochastic semiclassical gravity as well as methods to deal with them. 

Throughout the paper we use natural units with h = c = 1 and the (-f,-f,-)-) convention of Ref. ^|. We also 
make use of the abstract index notation of Ref. [3l|. Latin indices denote abstract indices, whereas Greek indices are 
employed when a particular coordinate system is considered. 

II. SEMICLASSICAL GRAVITY AND LINEAR RESPONSE THEORY 

A possible first step when addressing the interplay between gravity and quantum field theory is to consider the 
evolution of quantum matter fields (matter field is referred to here as any field other than the gravitational one) on 
a classical spacetime with a nontrivial geometry, characterized by a metric gab- As opposed to the situation for a 
Minkowski spacetime, there is in general no preferred vacuum state for the fields and particle creation effects naturally 
arise, such as Hawking radiation for black holes, cosmological particle creation and the generation of primordial 
inhomogeneities in inflationary cosmological models. Quantum field theory in curved spacetimes (QFTCST) is 
by now a well-established subject (at least for free fields and globally hyperbolic spacetimes) [s^lssf . 
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QFTCST is only an approximation in that the matter fields are treated as test fields evolving on a given spacetime. 
Einstein's theory requires that spacetime dynamics determines and is determined by the matter field. Thus one needs 
to consider the back reaction of the quantum matter fields on the dynamics of the spacetime geometry, which naturally 
leads to the semiclassical theory of gravity, where the evolution of the spacetime metric gab is determined by the 
semiclassical Einstein equation 



where gab is the spacetime metric, Ga6[g] is the Einstein tensor and the matter source corresponds to the renormal- 
ized expectation value of the stress tensor operator of the matter fields (a prime was used to distinguish it from 
that introduced below after absorbing some terms). Here, A is the renormalized cosmological constant, k ~ 8ttG, 
with G = i/rrip being the Newton constant and rrip the Planck mass; a and /3 are renormalized dimensionless cou- 
pling constants associated with tensors Aa^fg], _Bafc[(7] needed for the renormalization of the logarithmic divergences 
|8ll |. The expectation value of the stress tensor operator exhibits divergences which are local and state indepen- 
dent. Introducing a covariant regularization and renormalization procedure, those divergences can be absorbed into 
the cosmological constant, the Newton constant multiplying the Einstcin-Hilbert term and the gravitational action 
counterterms quadratic in the curvature. The finite contributions from those counterterms give rise to the covari- 
antly conserved tensors Aab and Bab which result from functionally differentiating with respect to the metric the 
terms / d^x^/^C°'^'^'^Cabcd and / d'^x^/—gE? respectively, where Cabcd is the Weyl tensor and R is the Ricci scalar. 
Those contributions were explicitly written on the left-hand side of Eq. (^) , but from now on will be included in the 
renormalized expectation value of the stress tensor operator so that the semiclassical Einstein equation becomes 



The field operators appearing in the stress tensor operator for the quantum matter fields are in the Heisenberg picture 
and satisfy the corresponding equation of motion, which coincides with the classical field equation for fields evolving 
on that spacetime. In particular, if we consider a free scalar field, the field operator in the Heisenberg picture will 
satisfy the corresponding Klein-Gordon equation for that geometry. 

Given a manifold M. and a metric gab which characterize a globally hyperbolic spacetime, and a density matrix 
p which specifies the state of the quantum matter fields on a particular Cauchy hypersurface, the triplet {M,gab,p) 
constitutes a solution of SCG if it is a self-consistent solution of both the semiclassical Einstein equation Q and the 
equations of motion for the quantum operators of the matter fields evolving on the spacetime manifold A4 with metric 
gab- Those operators enter in turn into the definition of the stress tensor operator appearing in the semiclassical 
Einstein equation. 

The semiclassical Einstein equation has been justified in at least two different ways. One possibility is to argue, by 
assuming a number of reasonable axioms, that it is the only consistent way to extend the classical Einstein equation 
to couple the quantum matter fields to a classical metric ,21. .33j . Alternatively, it can be derived by considering N 
free matter fields weakly coupled to the gravitational field in the sense that the gravitational coupling constant times 
the number of fields NG remains constant in the limit N ^ oo [3^ . The semiclassical Einstein equation can then be 
shown to correspond to the dynamical equation for the evolution of the expectation value of the metric in the limit 
of large N [s^l- Of course, in reality N is finite and the semiclassical Einstein equation can only be understood as the 
lowest order contribution in a expansion. 

Functional methods based on path integrals are useful not only in the rendition of ideas but also in actual computa- 
tions. However, the usual in-out formalism suitable for computing transition matrix elements in scattering problems 
is not appropriate to deal with back-reaction problems in which one is interested in the causal evolution of true 
expectation values from their initial values. The closed-time-path (CTP) forrnalism, which naturally yields true ex- 
pectation values and causal evolution equations for their dynamics 's^, 36, 37, 38, 39], should be used instead. This 
formalism has been applied to study a number of situations involving gravitational back-reaction effects of quantum 
fields |40l I41I fi^. I43L [4^ . (See also Ref. jl^ for an interesting application to an analogous situation in QED, where 
the back reaction of charged quantum fields on the dynamics of the expectation value of the electromagnetic field was 
considered.) This includes minisuperspace models which restrict the possible geometries to Robertson- Walker metrics 
and consider perturbative deviations from the massless conformal case for the matter fields E3, ^3 , as well as 
small metric perturbations (of a less restricted form) around self-consistent solutions of SCG [iy, llj . 

More specifically, given a background metric gab which is a solution of the semiclassical Einstein equation in SCG, 
one can compute the CTP effective action on the perturbed metric gab = gab + hab up to quadratic order in the 
metric perturbations hab- Taking the functional derivative of the CTP effective action with respect to hab, one gets 
the perturbed version of Eq. Q to linear order in the metric perturbations: 
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where the superindex (1) was used to denote that only terms hnear in hab should be considered. The linearized 
Einstein tensor G^^^j [g + h] conies from the Einstein-Hilbert term in the gravitational action. On the other hand, the 



contribution to the CTP effective action which results from functionally integrating the matter fields involving the 
expectation value and two-point correlation functions of the stress tensor operator on the background geometry yields 
the following result for the linearized expectation value of the stress tensor operator: 



-2(H 



2{M-hf'' (x), 



(4) 



where we have introduced the notation A - B = f d'^y^J—g{y)A°-^{y)Bab{y), and the kernels H and M are given by 



M 



abed I 



x,y) 



H^'^'^ix^y) = --lm[T*T'^'[^,g;x)T^^[^,g;y)] 

1 5((f-''[^,gab;x)] 



T^' [ip,g;x),T^'' [ip,g;y) 



Sgcdiy) 



(5) 
(6) 



where the notation T* was employed to indicate that the spacetime partial derivatives appearing in the time-ordered 
operators also act on the theta function implementing the time ordering. The functional derivative appearing on the 
right-hand side of Eq. (j^ should be understood to account only for the explicit dependence on the metric: the implicit 
dependence through the field operator ip[g] is excluded [s^. 

The previous result for the expectation value of the stress tensor when small metric perturbations around a back- 
ground solution of SCG are considered has been obtained in two ways: 1) by applying the usual techniques of linear 
response theory (see, for instance, j4^) to SCG j3, 47, 48J, and 2) by employing the influence functional formalism 
|49l ISOl to linearized metric perturbations regarded as an open quantum system interacting with an environment 
constituted by the quantum matter fields [5ll l52| . The influence functional approach also provides the noise kernel 
which underscores the stochastic nature of the dynamics for the metric perturbations. We will employ this method 
in the next Section. 

The explicit expression for the linearized expectation value (T^^^ '^''[g + h; x)) in the particular case of a Minkowski 
background spacetime and a free scalar field in the Minkowski vacuum state was obtained in Refs. and for 
a massless field and in Refs. for a field with an arbitrary mass. For our discussion we have included them 

in Appendix IbI where a global incrtial coordinate system {x^} for the Minkowski background is used. According to 

can be written entirely in terms of the linearized Einstein tensor G^^-*"^. Taking that into 
account, the expression for the linearized semiclassical Einstein equation (jJl in Fourier space becomes 



where 



with 



0) 



F2{p) = - y [Ha{p) + 3i?Bb) - 2a - 6^ 



(7) 
(8) 

(9) 
(10) 



where a (3 are some constants which include the renormalized parameters a and /3 in Eq. (^), and the kernels Ha{p) 
and H-b{p) are defined in Eqs. IjBSp and HB4|I of Appendix IbI The solutions of Eq. Q were analyzed in Refs. (l6| 
and |2l| for the massless case and Ref. for the general case. There is an obvious solution for G^^"* (p) — 0, which 
corresponds to the linear gravitational waves propagating in Minkowski spacetime. In addition, there are solutions 
of the form g\^2 oc 5{p^ ~ Pq) for particular values of (positive or negative) comparable to rn^. Since they exhibit 
characteristic timescales or lengthscales comparable to the Planck scale, where semiclassical gravity is not expected 
to be reliable anymore, they are usually regarded as unphysical. A more detailed discussion on this kind of solutions 
is given in Sec. II VI and Appendix IeI 

In fact, as will be explained in Sec. II VI if one quantizes the linearized metric perturbations, Eq. |7J coincides with the 
equation governing the evolution of the expectation value of the operator hab for the metric perturbations. Therefore, 
an analysis of the stability of the solutions of linearized semiclassical Einstein equation |(7J| can be equivalently 
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understood in terms of the evolution for the expectation value of hab- In Sec. II VI we will argue that a stability analysis 
for solutions of SCG with respect to small quantum corrections based solely on the expectation value of the metric 
perturbations is incomplete and should be extended to take into account the metric fluctuations as well. Stochastic 
semiclassical gravity is particularly well-suited to study the fluctuations of the metric and will play an important role 
in our later discussions. Therefore, in the next section we briefly review the formalism based on the Einstein-Langevin 
equation for small metric perturbations around semiclassical solutions within the framework of stochastic semiclassical 
gravity. 

III. EINSTEIN-LANGEVIN EQUATION FOR METRIC PERTURBATIONS AROUND A GIVEN 

BACKGROUND 

The semiclassical Einstein equation, which takes into account only the mean values, is inadequate whenever the 
fluctuations of the stress tensor operator are important. An improved treatment is provided by the Einstein-Langevin 
equation of stochastic gravity, which contains a (Gaussian) stochastic source with a vanishing expectation value 
and a correlation function characterized by the symmetrized two-point function of the stress tensor operator. This 
theory has been discussed by a number of authors P, |3, E3i 1^ EE US- Consider a globally hyperbolic 

background spacetime and an initial state for the quantum matter fields (one usually restricts to free fields) which is 
a self-consistent solution of SCG, i.e., it satisfies the semiclassical Einstein equation with the expectation value of the 
stress tensor operator obtained by considering the evolution of the matter fields on the same background geometry. 
The Einstein-Langevin equation governing the dynamics of the linearized perturbations hat around the background 
metric gab is given by 

Gi\^ [9 + h]^K (fil> [g + h]) +K U [g] , (11) 

\ / rcn 

where the Gaussian stochastic source ■Cah[5] is completely characterized by its correlation function in terms of the 
noise kernel Afabcd{x,y), which accounts for the fluctuations of the stress tensor operator, as follows: 

{^ab[g;x)^cd[g;y))^ =J^abcd{x,y) = ^{{iab[g;x),icd[g;y)}) , (12) 

where tab = Tab — (Tab) and (. . .) is the usual expectation value with respect to the quantum state of the matter 
fields, whereas (. . denotes taking the average with respect to all possible realizations of the stochastic source S^ab- 
Note that any local term quadratic in the curvature arising from finite contributions of the counterterms required 
to renormalize the bare expectation value of the stress tensor operator has been absorbed into its renormalized 
version {T^^l^g + /i])rcn- It should also be emphasized that solutions of the Einstein-Langevin equation for the metric 
perturbations are classical stochastic tensorial fields, not quantum operators. 

The precise meaning that should be given to these stochastic metric perturbations and the relation of the corre- 
sponding stochastic correlation functions to the quantum fluctuations resulting from quantizing these metric pertur- 
bations will be discussed below. Before doing so, it is, however, useful to mention some of the basic properties of 
the Einstein-Langevin equation (a more detailed discussion can be found in Refs. 0,0, |0|). First, when taking the 
average of the Einstein-Langevin equation with respect to all the possible realizations of the stochastic source, 
we recover the semiclassical Einstein equation Q, as follows straightforwardly from the vanishing expectation value 
of the stochastic source. Secondly, the integrability of the Einstein-Langevin equation is guaranteed, in the same way 
as in the semiclassical Einstein equation, by the conservation of the matter sources. The conservation of the expec- 
tation value of the stress tensor operator follows immediately from the fact that the divergence (with respect to the 
covariant derivative) of the stress tensor operator vanishes when the equation of motion of the matter field operators 
is satisfied {e.g. the Klein-Gordon equation for a scalar field). On the other hand, the fact that (V"Cab(a;))^ — 
and {£,ab{x)'^'^^cd{x))^ = 0, which completely characterize the Gaussian stochastic field W^^^abix), guarantees the 
conservation of the stochastic source (unless otherwise stated, from now on all the covariant derivatives are taken 
in the background metric and indices are raised and lowered using also the background metric). The previous two 
equalities are, respectively, a consequence of the vanishing expectation value of the stochastic source and the fact 
that V^tab{x) = 0. Finally, the Einstein-Langevin equation is invariant under gauge transformations corresponding 
to infinitesimal diffeomorphisms characterized by an arbitrary vector field Ci^), which generate transformations for 
the metric of the form hab hab + VaCb + '^bCa- This fact can be seen by realizing that the stochastic source does 
not depend on the metric perturbations, whereas the terms depending on the metric perturbations are all together 
gauge invariant. This is because they correspond to perturbing the semiclassical Einstein equation for the background 
metric ^ab, which is automatically satisfied since the background configuration under consideration for the metric and 
the state of the matter fields is a solution of SCG. 
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The Einstein-Langevin equation had been previously derived making use of a formal analogy with open quantum 
systems and employing the influence functional formalism |49ll5C| . This form was also justified in Ref. [s^ by arguing 
that it is the only consistent generalization of the semiclassical Einstein equation which takes into account the lowest 
order effects due to the fluctuations of the stress tensor operator. In fact, making use of a large N expansion, one 
can show that the stochastic correlation functions for the metric perturbations obtained from the Einstein-Langevin 
equation coincide with the leading order contribution to the quantum correlation functions in the large N limit. The 
details of the derivation will be given in Ref. and are summarized in Appendix |C1 for the particular case of a 
Minkowski background, to which we will restrict in the present discussion. In particular, the two-point stochastic 
correlation function is equivalent to the symmetrized quantum correlation function to leading order in provided 
that one also averages over the initial conditions for the solutions of the Einstein-Langevin equation distributed 
according to the Wigner functional characterizing the initial state of the metric perturbations (see Eq. |IC11|) in 
Appendix El for the definition of the Wigner functional). It is, therefore, convenient to express the solutions of the 
Einstein-Langevin equation as 

hab{x) = Si°)(a;) + A^(G,.et • Oab{x), (13) 

where R, = Nk is the rescaled gravitational coupling constant introduced in Appendix |0 '^^ab i^) ^ solution of the 
homogeneous part of the Einstein-Langevin equation (|ll|l containing all the information about the initial conditions 
(by homogeneous part we mean Eq. excluding the stochastic source, which coincides with the semiclassical 
Einstein equation 0), and Grct{x,x') is the retarded propagator with vanishing initial conditions associated with 
that equation (see Sec. IE 3l in Appendix IeI for important remarks on the propagator). Using Eq. H12|) . we can then 
get the following result for the symmetrized two-point quantum correlation function: 

\ {{Kb{x)Xa{x')}) = (Ei°)(x)Ef;(x'))^„ + % (Grot • A^- (G,ct)^),,,, (x,x'), (14) 

ab ' (i) 

where the Lorentz gauge condition V^{hab — l/2f?ah^o) = as well as some initial condition to fix completely the 
remaining gauge freedom of the initial state should be implicitly understood, and the stochastic source was rescaled 
according to Appendix lUl so that {(,ab[9', x)(,cd[g',y))^ = {^/N)J\fabcd{x,y), where J\fabcd{x,y) is the noise kernel for a 
single field. 

This result is analogous to that obtained in Ref. 59] for linear QBM models and briefly summarized in Appendix 1X1 
It should be emphasized that, similar to that case, there are two different contributions to the symmetrized quan- 
tum correlation function. The first one is connected to the quantum fluctuations of the initial state of the metric 
perturbations and we will refer to it as intrinsic fluctuations. The second contribution, proportional to the noise 
kernel, accounts for the fluctuations due to the interaction with the matter fields, and we will refer to it as induced 
fluctuations. In the next section we will formulate a generalized stability criterion for the solutions of SCG which 
involves the quantum fiuctuations of the metric. In particular we will see that the induced fiuctuations will play an 
important role on that issue. 

The noise kernel that we need for our discussions is for the particular case of a Minkowski background spacetime 
with a scalar field in the Minkowski vacuum. It was obtained in Ref. [s^l and is given by Eq. ljB7|) in Appendix IbI 

IV. STABILITY CRITERION FOR SOLUTIONS OF SEMICLASSICAL GRAVITY 

In this Section we will propose a criterion for analyzing the stability of a given solution of SCG with respect to 
small quantum corrections, associated with quantized metric perturbations around a background geometry. As an 
important example, we will apply this to the particular case of a Minkowski background with N scalar fields in the 
Minkowski vacuum state. 



A. Stability of Minkowski space: previous criteria 

The stability of metric perturbations around a Minkowski spacetime interacting with quantum matter fields in 
their Minkowski vacuum state was first studied in the context of SCG by Horowitz (16|. He considered massless 
conformally coupled scalar fields and found exponential instabilities for the linearized metric perturbations with 
characteristic timescales comparable to the Planck time. Those solutions are closely related to the higher derivative 
countertems required to renormalize the expectation value of the stress tensor operator (see, however. Appendix IeI for 
further comments on this point) and are analogous to the runaway solutions commonly present in radiation reaction 
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processes such as those considered in classical electrodynamics |60|, |61| . It is generally believed that the runaway 
solutions obtained by Horowitz are an unphysical artifact since they involve scales beyond the regime where SCG is 
expected to be reliable (in fact, this statement can be naturally formulated when regarding general relativity as a low 
energy effective theory). 

Since the existence of terms with higher derivatives in time implies an increase in the number of degrees of freedom 
(in an initial value formulation, not only the metric and its time derivative should be specified, but also its second 
and third order time derivatives), it seems plausible that, by restricting to an appropriate subspace of solutions of the 
semiclassical Einstein equation, one can reestablish the usual number of degrees of freedom in general relativity and, 
at the same time, get rid of all the unphysical runaway solutions. Following this line of thought Simon proposed that 
one should restrict to solutions which result from truncating to order h an analytic expansion in h (or equivalently in 
Ip, the Planck length squared) 19, 2Q|. Together with Parker he also introduced a prescription to reduce the order of 
the semiclassical Einstein equation which was computationally convenient in order to obtain solutions corresponding 
to such truncated perturbative expansions in h |22l | . 

On the other hand, Flanagan and Wald argued that Simon's criterion based on truncating to order h solutions 
which correspond to analytic expansions in h seemed too restrictive since it only allowed small deviations with respect 
to the classical solutions of the Einstein equations (21i] . In particular, one would miss those situations in which the 
small semiclassical corrections build up to give significant deviations at long times, such as those corresponding to 
the evaporation of a macroscopic black hole (with a mass much larger than the Planck mass) by emission of Hawking 
radiation. Furthermore, they illustrated with simple examples that there are cases in which one expects that no 
solutions of the semiclassical equation are analytic in h. Therefore, they suggested that, rather than trying to restrict 
the subspace of acceptable solutions, one should simply transform the semiclassical equation, by making use of Simon 
and Parker's order reduction prescription, to a second order equation which were equivalent to the original equation 
up to the order in h (or 1^) under consideration. All the solutions of the second order equation should then be regarded 
as acceptable, even if they are not analytic in h. Obviously, one could only extract physically reliable information 
from those solutions for scales much larger than the Planck length. 

Yet another prescription was proposed by Anderson, Molina-Pan's and Mottola 3] on the stability of small metric 
perturbations around the Minkowski spacetimc. They got rid of the unphysical runaway solutions by working in 
Fourier space and discarding those solutions which corresponded to 4-momenta with modulus comparable or larger in 
absolute value than the Planck mass. However, it is not clear how this procedure could be generalized to situations 
where working in Fourier space is not adequate, as in time-dependent background spacetimes. 

The consequences of both the order reduction prescription introduced by Simon and Parker and advocated by 
Flanagan and Wald, and the procedure employed by Anderson et at are rather drastic, at least when applied to the 
case of a Minkowski background, since one is just left with the solutions of the sourceless classical Einstein equation 
corresponding to linear gravitational waves propagating in Minkowski. In fact, the situation was not completely trivial 
for Flanagan and Wald, who were interested in analyzing whether the averaged null energy condition (ANEC) was 
satisfied in SCG by considering perturbations of the Minkowski solution, because they also perturbed the state of the 
matter fields. The order reduction prescription also seems to exclude those solutions which correspond to inflationary 
models driven entirely by the vacuum polarization of the quantum matter fields .62] , such as the trace anomaly driven 
inflationary model initially proposed by Starobinsky 63]. To keep this kind of models. Hawking, Hertog and Reall 
considered a less drastic alternative to deal with the runaway solutions l65j|. Their procedure, which is analogous 
to some methods previously employed in classical electrodynamics for radiation reaction problems 60], is based on 
discarding solutions which grow without bound at late times (see Appendix ^ for further discussions on this and 
related issues). 



B. Generalized stability criterion 

How does one characterize the quantum state of the metric perturbations? The first candidate is the expectation 
value for the operator associated with the perturbation of the metric, hat- In fact, using a large N expansion, 
Hartle and Horowitz showed that the semiclassical Einstein equation can be interpreted as the equation governing 
the evolution of the expectation value of the metric to leading order in [s^]- Taking that result into account, 
the study of the stability of a solution of SCG by linearizing the semiclassical Einstein equation with respect to small 
metric perturbations around that solution can be understood in the following way: Take an initial state for the metric 
perturbations with a small nonvanishing expectation value for the operator hab, let it evolve, and sec if the expectation 
value grows without bound. 

However, in addition to the expectation value of hab the state of the metric perturbations will also be characterized 
by its fluctuations. In fact, if there was no interaction with matter fields so that the state for the metric perturbations 
evolved unitarily, the set of quantum correlation functions (for the operator hab) evaluated at equal times would 
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completely characterize the quantum state of the metric perturbations '84| . Let us now suppose that the evolution of 
the expectation value is stable (i.e., that it does not grow unboundedly with time) or even that it vanishes for all times. 
It is clear that the semiclassical solution cannot be regarded as stable with respect to small quantum corrections if the 
fluctuations of the state for the metric perturbations grow without bound. Therefore, the stability criterion stated in 
Ref. ,3J should be generalized: one also needs to take into account the fluctuations. According to Ref. Q, a necessary 
condition for the stability of a solution of SCG requires that no gauge invariant scalar quantity constructed just 
from the linearized metric perturbation hat (which satisfles the semiclassical Einstein equation linearized around the 
semiclassical solution under consideration) and its derivatives grows without bound. This criterion can be interpreted 
as a condition on the stability of the expectation value of the operator hab for the state of the metric perturbations. 
We claim that, in addition, the n-point quantum correlation functions for the metric perturbations (starting with 
n = 2) should also be stable. Considerations based on gauge-invariant variables will not be necessary because we will 
be dealing with expressions where the gauge freedom has been completely fixed. 

As explained in Appendix lUl to leading order in the CTP generating functional for the metric perturbations 
exhibits a Gaussian form provided that a Gaussian initial state for the metric perturbations with vanishing expectation 
value is chosen. All the n-point quantum correlation functions can then be obtained, to leading order in from 
the two-point quantum correlation function. Furthermore, any of the two-point quantum correlation functions can in 
turn be expressed in terms of the symmetrized and antisymmetrized correlation functions (the expectation values of 
the commutator and anticommutator of the operator hab)- To leading order in the commutator is independent 
of the initial state of the metric perturbations and is given by 2iK{Grct{x' , x) — G,-et{x,x')). On the other hand, 
the expectation value of the anticommutator is given by Eq. (|14|) and is the sum of two separate contributions: the 
intrinsic and the induced fluctuations. 

The first contribution in Eq. H14|l to the correlation function for the metric perturbations involves the solutions of the 
homogeneous part of the Einstein-Langevin equation which actually coincides with the linearized semiclassical 
equation for the metric perturbations around the background geometry. Similarly, Grot corresponds to the retarded 
propagator (with vanishing initial conditions) associated with the linearized semiclassical equation. Thus, solving 
the perturbed semiclassical Einstein equation not only accounts for the evolution of the expectation value of the 
metric perturbations, which will exhibit a nontrivial dynamics as long as we choose an initial state with nonvanishing 
expectation value, but also provides nontrivial information, even for a state with a vanishing expectation value, about 
the commutator as well as the intrinsic fluctuations of the metric. This implies that the analysis about the stability of 
the solutions of SCG can also be used to determine the stability of the metric perturbations with respect to intrinsic 
fluctuations. 

The new observation we make here is that the induced fluctuations can be important as well. Both the retarded 
propagator and the solutions of the linearized semiclassical Einstein equation depend, through the kernel H, on the 
expectation value of the commutator of the stress tensor operator on the background geometry and on the imaginary 
part of its time-ordered two-point function. However, they do not involve the expectation value of the anticommutator, 
which drives the induced fluctuations. Furthermore, although the expectation value of the commutator and the 
anticommutator are related by a fluctuation-dissipation relation in some particular cases |5lL Is^ , that is not true in 
general and the induced fluctuations need to be explicitly analyzed. 

To sum up, when analyzing the stability of a solution of SCG with respect to small quantum corrections, one 
should also consider the behavior of both the intrinsic and induced fluctuations of the quantized metric perturbations. 
Whereas information on the stability of the intrinsic fluctuations can be retrieved from an analysis of the solutions of 
the perturbed semiclassical Einstein equation, the effect of the induced fluctuations is properly accounted for only in 
the stochastic semiclassical gravity framework based on the Einstein-Langevin equation. 

C. Stability of Minkowski space from our criterion 

We now turn to the application of the criterion proposed in the previous subsection to the particular yet important 
case of Minkowski spacetime. As explained there, the existing results in the literature can be interpreted as analysis 
of the stability of the expectation value of the operator associated with the metric perturbations (see, however, 
Refs. ESS El). On the other hand, we also need to include in our consideration the fluctuations, characterized by 
the two-point quantum correlation function. 

Before proceeding to analyze the two-point quantum correlation functions it is convenient to decompose the metric 
perturbations around Minkowski in the following way : 

hab = 4>Vab + (V(,Vfc) - TyafcQ) ^ + 2V(,W6) + hji^ , (15) 

where f " is a transverse vector and h"^^ is a transverse and traceless symmetric tensor, i.e., Vaf ° = 0, V^ft.^j'" = and 
{h'^'^)a = 0. Similarly, any vector fleld characterizing an inflnitesimal gauge transformation can be decomposed as 
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_ ya^ _j_ ya ^ where IS a transverse vector field. It is then clear that the vectorial and one of the scalar parts of 
the metric perturbation corresponding to and ip respectively, can be eliminated by choosing a gauge transformation 
such that = —v°- and C, = (this will also imply a change for 0: — > + 2nC). 

When the Lorentz gauge W^{hab — ^/^rjabh'^) = is imposed, we get the following conditions on the metric 
perturbations: Ov^ = and Wb(f> — (which implies </> = constant). Any vector field characterizing the remaining 
gauge transformations compatible with the Lorentz gauge satisfies the condition = 0, which implies DV^ = 
and V°nC = 0. We can see that a vectorial gauge transformation compatible with the Lorentz gauge can still be 
used to eliminate the vectorial part (now both w° and V"" must be solutions of the D'Alambertian equation). On 
the other hand, a scalar gauge transformation such that DC — -~4> — constant (this is always possible for Minkowski 
with a trivial -simply connected- topology) can be introduced to get (/> = 0. Moreover, an additional scalar gauge 
transformation compatible with the Lorentz gauge and leaving (j) invariant, which is characterized by a C which 
satisfies the D'Alambertian equation = (or, equivalently, C,{p) = for = p^p'^'q^j,^ ^ in Fourier space), can 
be used to eliminate those contributions to 'ijj which correspond to Fourier modes ipiv) with p^ = while leaving the 
remaining contributions unmodified. From now on we will assume that the Lorentz gauge has been imposed and that 
the additional gauge transformations just mentioned have been carried out so that we are left only with the tensorial 
components as well as those modes of the scalar component ij} with ^ in Fourier space. 

One could select the gauge mentioned in the previous paragraph imposing suitable conditions on the reduced Wigner 
functional characterizing the initial state for the metric perturbations; see Appendix IeI for some additional comments 
on this point. However, as explained in Appendix^ asymptotic initial conditions should be considered in order to 
get a finite result for the metric correlation functions. Therefore, rather than fixing the gauge for some initial state 
at some finite initial time, we will work in Fourier space implicitly assuming asymptotic initial conditions and fixing 
the gauge as described above. 

In order to analyze the two-point quantum correlation function for the metric perturbations, we will make use of the 
results mentioned in Sec. Illll and described in some more detail in Appendix lUl In particular, we will exploit the fact 
that the stochastic correlation functions obtained with the solutions of the Einstein-Langevin equation coincide with 
the quantum correlation functions for the metric perturbations. Moreover, according to Eq. H14|l . the symmetrized 
two-point quantum correlation function has two different contributions: the intrinsic and the induced fluctuations. 
We proceed now to analyze each contribution separately. 



The first term on the right-hand side of Eq. (|14|) corresponds to the fluctuations of the metric perturbations due to 
the fluctuations of their initial state and is given by 



where we recall that Y}^^ [x) is a solution of the homogeneous part of the Einstein-Langevin equation (once the Lorentz 
gauge has been imposed) with the appropriate initial conditions. 

As mentioned in Sec. IIIII and Appendix O the homogeneous part of the Einstein-Langevin equation actually 
coincides with the linearized semiclassical Einstein equation 0. Therefore, we can make use of the results derived in 
Refs. mil EH, which are briefly summarized in Appendix El As described there, in addition to the solutions with 
g']^^{x) = 0, there are other solutions that in Fourier space take the form c'^jllip) oc 5(p^ — p§) for some particular 
values of Pq, but they all exhibit exponential instabilities with Planckian characteristic timescales. 

In order to deal with those unstable solutions, one possibility is to employ the order reduction prescription. 
We are then left only with the solutions which satisfy g|jV' {p) = (see Appendix • The result for the met- 
ric perturbations in the gauge introduced above can be obtained by solving for the Einstein tensor in that gauge: 
^'abiP) ~ (l/2)p'^(^/^i/(p) — ^/'^Vp.i^^pip))- Those solutions for hfj^^{p) simply correspond to free linear gravitational 
waves propagating in Minkowski spacetime expressed in the transverse and traceless (TT) gauge. When substituting 
back into Eq. H16() and averaging over the initial conditions we simply get the symmetrized quantum correlation 
function for free gravitons in the TT gauge for the state given by the reduced Wigner function. As far as the intrinsic 
fluctuations are concerned, it is clear that the order reduction prescription is too drastic, at least in the case of 
Minkowski spacetime, since no effects due to the interaction with the quantum matter fields are left. The method 
employed in Ref. 's'l, although slightly different, yields the same result. 

A second possibility, proposed by Hawking et al. 64, 65], is to impose boundary conditions which discard the 
runaway solutions that grow unboundedly in time and correspond to a special prescription for the integration contour 
when Fourier transforming back to spacetime coordinates (see Appendix IeI for a more detailed discussion). Following 



1. Intrinsic fluctuations 




(16) 
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that procedure we get, for example, that for a massless conformally coupled scalar field with a = |85| and /3 > 
the intrinsic contribution to the symmetrized quantum correlation function coincides with that of free gravitons plus 
an extra contribution for the scalar part of the metric perturbations 4> which renders Minkowski spacetime stable but 
plays a crucial role in providing a graceful exit for inflationary models driven by the vacuum polarization of a large 
number of conformal fields (such a massive scalar field would not be in conflict with present observations because, for 
the range of parameters usually considered, the mass would be far too large to have observational consequences p^)- 



2. Induced fluctuations 

The second term on the right-hand side of Eq. (|14|l corresponds to the fluctuations of the metric perturbations 
induced by the fluctuations of the quantum matter flelds and is given by 

^ (Gret • . (G,et)^),,,, [x, x') - Nk^ (G,et ' U ' (G.-ct )^) ,,,, {x. x'), (17) 

where Afabcd{x, x') is the noise kernel accounting for the fluctuations of the stress tensor operator, and (Grct)afccrf(a;, x') 
is the retarded propagator with vanishing initial conditions associated with the integro-differential operator Labcd{x, x') 
deflned in Eq. (|C10|I of Appendix O 

As shown in Appendix the symmetrized two-point quantum correlation function coincides with the stochastic 
correlation function obtained from the solutions of the Einstein-Langevin equation. In fact, the contribution corre- 
sponding to the induced quantum fluctuations, given by Eq. (|17|) . is equivalent to the stochastic correlation function 
obtained by considering just the inhomogeneous part of the solution to the Einstein-Langevin equation: the second 
term on the right-hand side of Eq. (|13|l . Taking all that into account, it is clear that we can make use of the results 
for the metric correlations obtained in Ref. [H^ by solving the Einstein-Langevin equation (the homogeneous part of 
the solution was not considered there). In fact, one should simply take A'^ = 1 to transform our expressions to those 
of Ref. [H^l and, similarly, multiply the noise kernel in the expressions of that Reference by N so that they can be 
used here, which follows straightforwardly from the fact that we have N independent matter fields. 

The same kind of exponential instabilities in the runaway solutions of the homogeneous part of the Einstein-Langevin 
equation (the linearized semiclassical Einstein equation) also arise when computing the retarded propagator Grct- In 
order to deal with those instabilities, similar to the case of the intrinsic fluctuations, one possibility is to make use of 
the order reduction prescription. The Einstein-Langevin equation becomes then G^^"* = K^ab- The second possibility, 
following the proposal of Hawking et ai, is to impose boundary conditions which discard the exponentially growing 
solutions and translate into a special choice of the integration contour when Fourier transforming back to spacetime 
coordinates the expression for the propagator. In fact, it turns out that the propagator which results from adopting 
that prescription coincides with the propagator that was employed in Ref. |5J|. However, it should be emphasized 
that this propagator is no longer the retarded one since it exhibits causality violations at Planckian scales. A more 
detailed discussion on all these points can be found in Appendix IeI 

Following Ref. [s^ , the Einstein-Langevin equation can be entirely written in terms of the linearized Einstein tensor 
G^^} (p) as follows: 

F^.aMP)G(^^"^(p) -^ep.(p), (18) 

which simply corresponds to adding the stochastic source to the linearized semiclassical Einstein equation Q, where 

Ffj,vai3{p) was given by Eq. One can then solve the stochastic equation for &,}}{p) and obtain its correlation 
function [s^ : 

(GW(p)G«(y))e = -n^D,.M{C^{p)V'{j>'))iD,.,8{p') 

= ^D^.MU^^'^'{j>)Dpa^,s{-p){2^f5{p + P% (19) 

In the last equality we have taken into account translational invariance. The noise kernel ^["^^^{p) is given by 
Eq. (|lj7|l in Appendix IBl and Dfj^^a/iip) is the propagator that results from inverting F^^ajiip) (see Appendix^ for a 
discussion on the uniqueness of this propagator) and is given by 

D^^uM - -^^V,(.Vp). - F,(p)f3(p) P'^^-^"^' (20) 
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with = rjfj^i, — Pp.Pu/p'^, Flip) ^'^d P2{p) given by Eqs. ® and and F3(p) = Fi{p) + Sp'^F2{p). On the other 
hand, if we make use of the order reduction prescription, we get 

(GW(p)G«b'))e = -^'{C'(.p)C'ip'))i = ^J^'^f'''\p)i2n)*S{p + p'). (21) 

Note that G^i^(p) is gauge invariant when perturbing a Minkowski background because the background tensor G^f^ 
vanishes and, hence, jC^Gf^j also vanishes for any vector field (. 

Finally, using the expression for the linearized Einstein tensor in the Lorentz gauge, G^^j] — {l/2)p'^h^i, with 
(l/2)77p,,/i", we obtain the correlation function for the metric perturbations in that gauge: 

{hAp)hpAp'))i = ^j^D^.MpW''^''\p)Dp.ys{^p)i27:rS{p + p'), (22) 

or 

{hAp)KAp'))i = ^j^Af^.pAp){^^)'S{p + p'), (23) 

if the order reduction prescription is employed. It should be emphasized that, contrary to the linearized Einstein 
tensor G^^\ the metric perturbation hab is not gauge invariant. This should not pose a major problem provided that 
the gauge has been completely fixed. 

The correlation functions in spacetime coordinates can be easily derived by Fourier transforming Eqs. (|22|l or H23|l . 
However, there is apparently an infrared divergence at p'^ = 0, at least for the massless case. For the massive case 
the result is finite because the noise kernel J\f"^^^{p) is proportional to 6'(— — Aw?), so that m? > guarantees 
that p^ — lies outside the domain of integration. On the other hand, in the massless case the terms of the form 
PpPuPpPa I [p^)"^ appearing when substituting the noise kernel in Eqs. (|22l) and (|23|l . give rise to infrared divergences 
when computing the Fourier transform. In fact, even if we exclude the massless case, the result would be finite, but 
it would become larger and larger as we chose a positive but arbitrarily small mass. In any case, such an infrared 
divergence seems to be just a gauge artifact [SG^ . 

We can conclude that, once the instabilities giving rise to the unphysical runaway solutions have been properly dealt 
with, the fluctuations of the metric perturbations around the Minkowski spacetime induced by the interaction with 
quantum scalar fields are indeed stable (if instabilities had been present, they would have led to a divergent result 
when Fourier transforming back to spacetime coordinates). It should be emphasized that no ultraviolet divergences 
related to the coincidence limit of the noise kernel appeared in the previous analysis because we implicitly assumed 
asymptotic initial conditions when working in Fourier space, as explained in Appendix^ Furthermore, in contrast 
to the intrinsic fluctuations, even when using the order reduction prescription there is still a nontrivial contribution 
to the induced fluctuations due to the quantum matter fields. 



V. DISCUSSION 



In this paper we make the point that an analysis of the stability of any solution of SCO with respect to small 
quantum corrections should consider not only the evolution of the expectation value of the metric perturbations 
around that solution, but also their fluctuations, encoded in the quantum correlation functions. Making use of a large 
N expansion, where N is the number of matter fields, the symmetrized two-point quantum correlation function for 
the metric perturbations can be decomposed into two distinct parts: the intrinsic fluctuations due to the fluctuations 
of the initial state of the metric perturbations itself, and the fluctuations induced by their interaction with the matter 
fields. The stability of the first contribution turns out to be closely related to the stability of linearized perturbations 
of the semiclassical Einstein equation, whereas the second contribution is equivalent to the stochastic correlation 
functions in stochastic semiclassical gravity obtained from solutions of the Einstein-Langevin equation. 

As a specific example, we analyzed the two-point quantum correlation function for the metric perturbations around 
the Minkowski spacetime interacting with N scalar fields initially in the Minkowski vacuum state. Once the ultraviolet 
instabilities (discussed in Appendix IE|) which are ubiquitous in SCG and are commonly regarded as unp hysical, have 
been properly dealt with by using the order reduction prescription or the procedure proposed in Refs. [6J,|63 , both the 
intrinsic and the induced contributions to the quantum correlation function for the metric perturbations are found to 
be stable. In fact, one gets an infrared divergence for the massless case when computing the inverse Fourier transform 
for the induced contribution to the correlation function of the metric, but that seems to be purely a gauge effect, as 
argued in footnote [86]. 
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The symmetrized quantum correlation function for the metric perturbations obtained is in agreement with the real 
part of the propagator obtained by Tomboulis in Ref. using a large N expansion (sj (he actually considered 
fermionic rather than scalar fields, but that just amounts to a change in one coefficient). Tomboulis used the in-out 
formalism rather than the CTP formalism employed in this paper. Nevertheless, his propagator is equivalent to the 
time-ordered CTP propagator when asymptotic initial conditions are considered because in Minkowski spacetime 
there is no real particle creation and the in and out vacua are equivalent (up to some phase which is absorbed in 
the usual normalization of the in-out propagator). The use of a CTP formulation is, however, crucial to obtaining 
true correlation functions rather than transition matrix elements in dynamical (nonstationary) situations (such as in 
an expanding Robertson- Walker background geometry), where the in-out scattering matrix might not even be well 
defined at all. 

As we pointed out in the Introduction, Ford and collaborators have stressed the importance of the metric fluctuations 
and inves tiga ted some of their physical iinplications 0, IE IE S IS S ^i 0| • They have considered both intrinsic 
[E IE ^1 ^1 and induced fluctuations 0, IE IE S 1 which they usually refer to as active and passive fluctuations 
respectively. However, they usually consider these two kinds of fluctuations separately and have not provided a 
unified treatment where both of them can be understood as different contributions to the full quantum correlation 
function. Moreover, they always neglect the nonlocal term which encodes the averaged back reaction on the metric 
perturbations due to the modified dynamics of the matter fields generated by the metric perturbations themselves 
|88|. Their justification is by arguing that those terms would be of higher order in a perturbative expansion. That 
is indeed the case when considering a Minkowski background if the order reduction prescription is employed, but it 
is not clear whether it remains true under more general conditions. In fact, as mentioned in Ref. for the usual 
cosmological inflationary models the contribution of the nonlocal terms can be comparable or even larger than that 
of the remaining terms. Finally, in order to deal with the singular coincidence limit of the noise kernel, in Ref. 
Ford and collaborators opted to subtract a number of terms including the fluctuations for the Minkowski vacuum. 
Even when no such subtraction was performed (because a method based on multiple integrations by parts was used 
instead) [EEEIlS' they usually discard the fluctuations for the Minkowski vacuum. Therefore, the information on 
the metric fluctuations around a Minkowski background when the matter fields are in the vacuum state is missing in 
their work. 

We close this Section by recalling a couple of partially open issues for which either a better understanding or a 
better treatment would be desirable. The first issue is the singular coincidence limit for the noise kernel. It seems 
clear that, when properly treating the noise kernel as a distribution, a finite result for the metric correlation function 
is obtained except for some divergent boundary terms at the initial time. There is a natural physical interpretation: 
the completely uncorrelated initial state that was considered becomes pathological when the number of modes of the 
environment is infinite. A simple way to overcome this problem and obtain a finite result for the correlation function is 
to switch on the interaction smoothly so that the modes of the environment with arbitrarily high frequencies become 
correlated with the system in a nonsingular way. However, in order to preserve the conservation of the source in 
the Einstein-Langevin equation, which guarantees the integrability of the equation through the Bianchi identity, the 
interaction has to be turned on adiabatically and asymptotically past initial conditions are required. Therefore, other 
procedures should be devised to address situations that require specifying the initial conditions at a finite initial time. 

The other question which deserves further study is the procedure employed to deal with the runaway solutions 
discussed in Appendix IeI The order reduction prescription is rather drastic as its net outcome is to discard entirely 
the contribution from the dissipation kernel (as far as the expectation value and the symmetrized two-point correlation 
function are concerned), which encodes the averaged back reaction of the matter fields on the metric perturbations. 
As for the method employed by Hawking et al. in Refs. [63. l65l |. we find the fact that the choice of the physical 
solutions at a given instant of time depends on the far future somewhat unsatisfactory, and discarding solutions which 
grow unboundedly in time could get rid of other possible instabilities which are physically meaningful. Furthermore, 
it is not clear whether both procedures could be implemented in a general case. 

To gain insight into some of the previous aspects, an interesting possibility is to consider an analogous situation in 
QED with the electromagnetic field regarded as an open quantum system interacting with an environment constituted 
by the charged quantum fields. In fact, the analogy between SCG and the equation for the expectation value of the 
electromagnetic field to leading order (order 1) in a large N expansion for N charge d q uantum fields has been 
discussed by a number of authors \^ 0, |3 . One step further was considered in Ref. (45j, where the evolution of 
the expectation value of the electromagnetic field was considered to next to leading order in 1/N (order 1/N). The 
two-point quantum correlation functions (the CTP propagators) for the electromagnetic field to leading order in 1/iV 
(order 1/N), which play a crucial role there, are completely analogous to the quantum correlation functions for the 
metric perturbations considered here. 
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APPENDIX A: INTRINSIC AND INDUCED FLUCTUATIONS IN A SIMPLE QUANTUM BROWNIAN 

MOTION MODEL 

In this Appendix we illustrate the importance of the fluctuations induced by the environment when considering the 
quantum fluctuations for an open system. As an example we will use a simple model which was analyzed in some 
detail in Ref. [soj : a linear quantum Brownian motion (QBM) model that consists of a harmonic oscillator, which will 
be referred to as the system, bilinearly coupled to a set of harmonic oscillators, which constitute the environment. 

In Ref. '5^ it was shown that a stochastic description based on a Langevin type equation could be used to gain 
information on the quantum properties of the open system. In particular, the symmetrized two-point quantum 
correlation function for the system turns out to be equivalent to the correlation function obtained in the context of 
the stochastic description: 

\ {{x{t,),x{h)}) = [{X{t,)X{t^))^)^^ , (Al) 

where (...) denotes the expectation value with respect to the quantum state of the system, x{t) is the position 
operator for the system in the Heisenberg picture, (. . .)^ denotes the average over all possible realizations of the 
stochastic source ^(t) and (. . ■) Xi,pi is the average over all possible initial conditions for the solutions of the Langevin 
equation distributed according to the reduced Wigner function for the initial state of the system. The functions 
X(t) appearing inside the stochastic averages are solutions of the Langevin equation L • X = ^, where L{t,t') — 
M{cP /dfi + rij^(,„)(5(i — t') + Hicn{t,t'), with i?icn being the renormalized kernel appearing in the real part of the 
influence action and • = J^*^ dt throughout this Appendix. Here ^(i) is a Gaussian stochastic source with vanishing 
expectation value and correlation function {£,it)£,{t'))^ = N{t,t') where N{t,t') is the noise kernel, being the kernel 
appearing in the imaginary part of the influence action (see Ref. |59l | for further details). When the environment is 
initially in a thermal equilibrium state, the noise kernel is explicitly given by N{t, t') = dioI{uj) coth /3cij cosuj{t—t'), 
where I{oj) is the spectral density function, which characterizes the frequency distribution of the oscillators in the 
environment. 

The solution of the Langevin equation can be written as X{t) — Ao(i) -I- (Grot • C){t): where Ao(t) is a solution of 
the homogeneous part of the Langevin equation which contains all the information about the initial conditions and 
Gret(i,^') is the retarded propagator with vanishing initial conditions. Substituting the previous expression for X(t) 
into Eq. HA1|I and using the properties of the stochastic source, one obtains the following result for the two-point 
quantum correlation function: 

i {{X[h),x{h)}) - (Ao(ii)Ao(t2))x^,p. + (Crct -A^- (Gret)^) {hM). (A2) 

where the first contribution corresponds to the intrinsic fluctuations connected to the dispersion of the initial state of 
the system, and the term proportional to the noise kernel reflects the fluctuations induced by the system's interaction 
with the environment. Note the close analogy between Eq. (|A2p and the expression for the symmetrized two-point 
quantum correlation function in the gravitational case, given by Eq. (|14|) . 

Let us specialize to the case of an ohmic environment, i.e., the case in which the spectral distribution function for 
the frequencies of the oscillators in the environment is of the form I{uj) = M^us, where M is the mass of the system 
harmonic oscillator and 7 is some constant proportional to the square of the system-environment coupling constant. 
Then the kernel i?ion becomes local with Hrcn{t, t') — M^5'{t — t') and the homogeneous solution Ao(t) takes on the 
following simple form: 



Ao(t) = e 



X.i cos Cl{t - tj) + ( 4- ^X, ) sin n{t - U) 



(A3) 



where Q, = -y/ri^^j^ — (7/2)^ and we considered the underdamped case (ilrcn > 7/2). A similar result also holds for the 
overdamped case with the trigonometric functions replaced by the hyperbolic functions. Due to the exponential factor, 
Ao(t) and hence the intrinsic fluctuations will decay at times much larger than the relaxation time 27"^ jsOj . In fact. 
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if we take the limit ti — > —oo, the contribution to the two-point correlation function from the intrinsic fluctuations 
completely vanishes and one is just left with the induced fluctuations. If the initial state of the environment was 
a thermal state, the dissipation kernel (the antisymmetric part of Hycn{t,t')) and the noise kernel are related by a 
fluctuation-dissipation relation which characterizes the balance between the noise induced by the environment and 
the dissipation effect so that the two-point correlation function remains bounded in time. 

From the example employed in this Appendix, it is clear that the induced fluctuations play an important role when 
considering correlation functions in open quantum systems. In fact, for asymptotically past initial conditions they 
become the entire contribution to the correlation function since the intrinsic fluctuations are completely damped by 
the dissipation. In a more general context, such as the gravitational case, the dissipation kernel will not damp the 
intrinsic fluctuations, but the induced ones will still play an important role. 



APPENDIX B: DISSIPATION AND NOISE KERNELS IN MINKOWSKI SPACETIME 

In this Appendix we provide the expressions for the dissipation and noise kernels of a free real scalar field when 
a Minkowski background spacetime is considered and the state of the fields is the Minkowski vacuum. The details 
of their derivation can be found in Ref . 54] . All the expressions in this Appendix are given in Fourier space and 
derived by making use of the translational invariance in terms of the inertial coordinates employed for the Minkowski 
background. Given any expression A{p)^ the corresponding expression in spacetime coordinates A[x — y) can be 
simply obtained from 



A{x - y) = 



(27r)4 



■■Aip). 



(Bl) 



The linear combination —2{M + H) with M and H given by Eqs. ©-® is commonly referred to as the polarization 
tensor in analysis based on linear response theory. We will use the term "dissipation kernel" , by analogy with the 
usual terminology employed in the context of open quantum systems jooj . In the case under consideration the local 
kernel Mlf^'^^{x — y) is proportional to the Einstein tensor and can be absorbed in a finite renormalization of the 
gravitational coupling constant. The expression in Fourier space for the nonlocal kernel HI^J^^^{x — y) is 



(B2) 



where a and (3 are constants which include the renormalized parameters a and (3 appearing in Eq. P^^i, is the 
projector orthogonal to p^, given by P^^, = rj^i, — p^^p^/p^ , and 



Ha{p) = 



Hb{p) = 




-Z7rsignp° 9{~p^ - Am^) \ 1 + 4— + ip{p^) 



-i7rsignp° 9{~p'^ - Arn?) \ 1 + A^ + if{p^ 



P 



3 p2 



1 m"^ 



(B3) 



where ^ is the parameter characterizing the coupling of the scalar field to the spacetime curvature through a term of 
the form —{^/2)R(fP in the matter Lagrangian, and (p{p'^) is given by 



tp{p ) = da In 
Jo 



1 + Aja(l-a) 




1 



'2\ —1 — 4^5" arccotan 

V p 



1 



p 



m 

-1-4— 



(B5) 



Using the renormalized version of Eq. Q in Fourier space, the dependence on the metric of the renormalized expec- 
tation value of the stress tensor operator can be written entirely in terms of the linearized Einstein tensor as follows: 



f^-'[g + h;p) 



= 2P^'^ ( -^HAip) + la- Hb{p) + 2p\ (g^^))" {p) + {Ha{p) - 2a) (g^^))''' . (B6) 
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Following Ref. |5J|, we have employed a renormalization scheme in which the renormalization scale is fixed to 
fi^ — m? . This is, of course, not possible for the massless case. Nevertheless, the expression for the massless 
case can still be obtained by adding a term h\{w? / ^'^) to Eq. (|B5p and subtracting (19207r^)~^ ln(/x^/TO^) and 
(5 — l/6)^(967r^)~^ ln(^^/m^) respectively from 2a and 2/3 in Eqs. (|B2p and (|B6p before taking the hmit 0. 
The renormalized parameters will then depend on the arbitrary scale /i. If desired, it is always possible to choose 
a = by fixing the renormalization scale to some appropriate value. 
Finally, the expression for the noise kernel in Fourier space is given by 




(B7) 



APPENDIX C: STOCHASTIC AND QUANTUM CORRELATION FUNCTIONS 

It was initially believed that some kind of environment-induced decoherence mechanism was required to realize 
the stochastic dynamics described by the Einstein-Langevin equation [sills^ . Later, in Ref. j5^] it was shown that, 
even in the absence of decoherence, a stochastic description based on a Langevin type equation contains nontrivial 
information on fully quantum properties of simple linear open quantum systems. In particular, the reduced Wigner 
function of the system (see, for instance, Ref. jZli] for the definition and properties of the Wigner function) can be 
expressed as a double average for the solutions of the Langevin equation with respect to both the different realizations 
of the stochastic source and the initial conditions, which are distributed according to the reduced Wigner function 
at the initial time. This expression can then be used to derive the master equation governing the time evolution of 
the reduced Wigner function (or, equivalently, the reduced density matrix). Furthermore, the stochastic correlation 
functions for the solutions of the Langevin equation are actually equivalent to quantum correlation functions for the 
system observables. 

Although the previous results were obtained in Ref. 59] for linear open quantum systems, they can be extended to 
the case of nonlinear quantum field theories provided that some kind of Gaussian approximation for the corresponding 
influence functional is considered. In fact, in Ref. [s^ it will be explained in detail how those results can indeed be 
shown to hold for the metric fluctuations around a given background spacetime by properly treating the gauge freedom 
and the corresponding dynamical constraints. More precisely, when considering N free quantum matter flelds weakly 
interacting with the gravitational field in the sense that the gravitational coupling constant times the number of flelds 
remains constant in the large N limit, the stochastic correlation functions can be shown to coincide with the leading 
order contribution to the quantum correlation functions of the metric perturbations in a large iV expansion. 

Here we briefly sketch, in the context of a large expansion, some of the key aspects in the derivation of the result 
stated above. The details will appear in Ref. [Sg and were partially included in Ref. j72i |. 

We will consider metric perturbations around a globally hyperbolic background spacetime (that will be specialized 
to Minkowski spacetime at some point) regarded as an open quantum system interacting with the quantum matter 
flelds, which constitute the environment. In particular we will consider N minimally coupled free scalar flelds, but 
the main result can be generalized to nonminimally coupled scalar flelds or even vectorial and fermionic fields. The 
action for the combined system is the sum of the gravitational action Sg plus the action for the matter fields Si^. The 
gravitational action is given by the usual Einstein-Hilbert term, the corresponding boundary term (which should be 
included to have a well-defined variational problem and will later be important) and the usual counterterms required 
to renormalize the divergences arising when functionally integrating the matter fields: 

Sg^— d'^xy/^R{g) + — d^xVffs^i^a (5) + (counterterms), (CI) 

2k Jj^ k Js=dM 

where gab = 9ab + hab is the perturbed metric, gab is the background metric and the gravitational coupling constant 
K — STr/rrip was rescaled to R/N so that the product of the rescaled gravitational constant times the number of fields 
remains constant in the limit iV — > 00. The action for the matter fields is 

N ^ 

= - I] / d^^V^- {r'^a^j^b^j + mVj) , (C2) 

where m is the mass of the scalar field. In fact, we will not take the limit N 00, but rather use the expansion in 
1/A'^ as a useful way to organize our computation and the contributions that are included. At the very end one can 
always substitute back the rescaled gravitational constant in terms of the physical one. 
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The closed time path (CTP) generating functional for the metric perturbations, from which true expectation values 
and correlation functions can be obtained [s^, [s^, IH, ^3, 0| , is given by 

N 

where we have used the notation A - B = J d^yy^—g{y)A°-''{y)Babiy)- The density matrices for the initial state of the 

fields and the metric perturbations, which are all assumed to be initially uncorrelated, are p[(p'j'\ Vj^*"*] ^^'^ Pi^ab^ ^'Jd ^] 
respectively. The gauge freedom in the path integrals for the metric perturbations should be properly treated, as 
briefly described below. 

The first step is to integrate out the matter fields using the influence functional formalism of Feynman and Vernon 
for open quantum systems |49l l50| . The influence action Sif is defined as 



N 

jsMh,h'] _ 



Up to quadratic order in the metric perturbations it is given by [HH IH^ I 

5iF [Safe, Aafc] = iV f Z • A + A • (if + A/) • S + • • aV (C5) 



where we have introduced the semisum and difference variables Yiab — {hab + h'^^b) ^^^d Aab — h'^b ~ ^ab, Z"-^{x) = 
-(l/2)(T°^[(^,g;x)) and the kernels H, M and A/" were defined in Eqs. ©, ® and As explained in SeclUl the 
kernels H and M exhibit divergences that are canceled by renormalizing the gravitational coupling constant and the 
cosmological constant in the bare gravitational action as well as the coupling constants of the counterterms quadratic 
in the curvature. We will not need terms of higher order in the metric perturbations because they give contributions 
to the connected part of the CTP generating functional (given by W^^^ = —i\nZ^^^) of higher order in 1/A^. 
This is also true for the terms in the gravitational action 5g, which implies that we do not have to consider graviton 
vertices. In order to show that, when computing the connected part of the CTP generating functional to leading order 
in 1/A^, it is indeed sufficient to keep just those terms in the gravitational action and the influence action which are 
at most quadratic in the metric perturbations, one can first compute the generating functional resulting from that 
approximation, Z^r^ , and then show that including terms of higher order in the metric perturbations would yield 
corrections to the connected part of the generating functional of higher order in 1/A^. 

We now specialize to the case of a Minkowski background, i.e., gab — Vab, and consider a family of Cauchy hypersur- 
faces which foliate the spacetime into constant time hypersurfaces of a given incrtial frame in Minkowski spacetime. 
The initial states for the metric perturbations and the matter fields are specified on one of these hypersurfaces, which 
will be denoted by Si. Another hypersurface is chosen as the final hypersurface Sf so that any spacetime region 
of interest lies between them. Then we integrate by parts those contributions to the Einstein-Hilbert term of the 
gravitational action involving two derivatives acting on the same factor and impose the Lorentz gauge condition 
V^/iofc — (we recall that the indices are raised and lowered using the background Minkowski metric, all the covariant 
derivatives are taken in this background metric and Aab = Aab — (l/2)'7a6^a)- The boundary terms resulting from 
integration by parts are canceled by the boundary terms included in the gravitational action 91] and the expression 
Sg[h, h'] ~ Sg[h] — Sg[h'] up to quadratic order in the metric perturbations becomes 

5g[Eafc,A,d] = — / d^^y^VaA^^V^Efcc + (counterterms). (C6) 

Next, we introduce the momentum canonically conjugate to Aab |9^ . which is given by 

where we employed the notation Aab = n'^'^ c^ab for the covariant derivative with respect to the normalized and 
future-directed timelike vector orthogonal to the family of Cauchy hypersurfaces including the initial and final 
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hypersurfaces Si and Sf. Finally, one can integrate again by parts so that 

4^ Jm JSf=dM JS,=dM 

+ (counterterms), (C8) 

where the indices (i) and (/) denote quantities evaluated on Si and 5/ respectively. It should be emphasized that 

^ib ^\ib simply correspond to the spacetime metric evaluated on those hypersurfaces and should not be confused 
with the induced metric in the usual ADM formulation. Furthermore, the contribution for Sf will not be relevant 
because, when computing the CTP generating functional, we should take hab ~ h'^i, (which implies Aab — 0) on the 
final hypersurface. 

Changing to the new current variables Jj"^ — {Jab + J'ab)/"^ ^"^^ •^ab ^ ^afc ^ '^ab, and functionally integrating with 
respect to Acd^ one gets the following expression for the generating functional: 

where the functional integral with respect to Safe is restricted to those configurations that satisfy the Lorcntz gauge 
condition, and we introduced the integro-differential operator 

Labcd{.x,x') = (l/2)(ry,er/6d - r]abr]cdmU5{x ~ x') + 2fiHiZ''J {x - x') + 2kM^,Z'{x - x'). (CIO) 

K is some normalization constant which can be eventually determined by demanding that ^'tp I^oi' "^cdl ~ ^ when 
we take J^f, = — 0, and Wr[s|(*^\ is the reduced Wigner functional for the metric perturbations at the initial 
time, which is defined in terms of the reduced density matrix at the initial time as 

T^,.[E«,n^,f)] = (1/2^) j dA«exp(ffl^,^A«)p,[S«,A«]. (CU) 

The Lorentz gauge condition does not fix completely the gauge freedom under local diffeomorphisms. However, it can 
be completely fixed by imposing additional gauge fixing conditions on the state of the metric perturbations at the 
initial time. Thus, from now on it should be understood that some appropriate condition such as the transverse and 
traceless gauge has been imposed on the reduced density matrix (or, equivalently, the reduced Wigner functional) at 
the initial time. 

Introducing a suitable functional change, the CTP generating functional can be rewritten in the following way when 
taking Jfb = 0: 

4TP^[jfb = 0,J.^J = ((e'^^-^) \ , (C12) 

ah ' (i) 

where the expectation values (. . and (. . are defined as 

= / dsSdn^f)...T4^r[s«,n^f)], (C13) 

(. . .)^ = (27rAA/7V)-i/2 j VU ■ ■ ■ e-^«-^"-«. (C14) 

The 12ab{x) inside the expectation values in Eq. ljC12p satisfies the equation 

{L ■ J:)ab{x) = kU{x), (C15) 

with initial conditions E^'^^ and = — (4K/A^)n(f) on the initial hypersurface Si. From Eq. (|C14|I it becomes clear 
that one can formally interpret ^ab as a Gaussian stochastic source with a vanishing expectation value and whose 
correlation function is given by the noise kernel. Eq. IjClSp can then be regarded as a stochastic Langevin equation 
and coincides with the Einstein-Langevin equation expressed in the Lorentz gauge (when integrated with the metric 
perturbation hab, the first term on the right-hand side of Eq. (jClOp corresponds to the linearized Einstein tensor, 

whereas the last two terms correspond to {Tj^l\g + h]), as follows from Eq. Q). Furthermore, 4tp ["^ai ~ '^cdl 
is also the generating functional for the stochastic correlation functions for the solutions of the Einstein-Langevin 
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equation and, therefore, the stochastic correlation functions are actually equivalent to quantum correlation functions 
for the metric perturbations. 

The solutions of Eq. (|CT5)l can be expressed as 



(C16) 



where S^^^'' is a solution of the homogeneous equation {L ■ T,)ab = 0, which coincides with the linearized semiclassical 
Einstein equation for the metric perturbations in the Lorentz gauge, with all the information on the initial conditions, 
and {G ret) abcdix, x') is the retarded propagator associated with the integro-differential operator Latcdix, x') with 
vanishing initial conditions on Si. The CTP generating functional for a nonvanishing J^j, can then be written as 
follows: 



(C17) 



where we introduced the notation {A^)abcd{x,y) = ^cda&(y,a;). 

It is interesting to consider the particular case of the two-point correlation function. Functionally differentiating 
twice with respect to and then taking and J^j^ equal to zero, one gets the following result for the symmetrized 
quantum correlation function for the metric perturbations: 



^hab{x),hcd{x')^ 



+ -(G,efAA.(G,etn^,^,(x,x') 



(C18) 



One can see that there are two separate contributions to the two-point correlation function: the first one is related to 
the dispersion of the initial state for the metric perturbations, whereas the second one is proportional to the noise kernel 
and accounts for the fluctuations induced by their interaction with the environment (in this case, the quantum matter 
fields). We refer to these two contributions as intrinsic and induced fluctuations respectively. Furthermore, taking into 
account Eq. (|C12p . we see that, under the aforementioned conditions, the symmetrized quantum correlation function 
for the metric perturbations is equivalent to the stochastic correlation function obtained in stochastic semiclassical 
gravity by solving the Einstein-Langevin equation. 

From the expression for the generating functional in Eq. ljC17ll one can get the remaining two-point quan- 
tum correlation functions to leading order in 1/N. In particular the commutator is given by {[hab{x),hcd(x')]) = 
2iK{G-cet{x' , x) — Grct{x, x')), and by combining the commutator and the anticommutator the rest of two-point func- 
tions can be easily obtained. Moreover, assuming a Gaussian initial state with vanishing expectation value for the 
metric perturbations, the expression for the generating functional in Eq. ljC17p becomes Gaussian and any other 
n-point quantum correlation function has a simple expression in terms of the two-point functions. 

The exact CTP generating functional is given by 



ZcTp[Jab, Jed] = exp iS; 



7(LO) 
^CTP 



,6 J '^cdlj 



(C19) 



where S'int Pab, ^cd] corresponds to all the terms in the gravitational action or the exact influence action of cubic or 
higher order in the metric perturbations. In order to consider and evaluate the different contributions to Eq. ljC19|l . it 
is convenient to introduce the corresponding Feynman rules and diagrams (in the CTP formulation) as follows: each 
term in S'int [^/^J^b' ^/^Jcd] gives rise to a vertex with the same number of legs as the total power of the functional 
derivatives S/6J^f^ and S/6J^ appearing in that term, and the CTP propagators simply correspond to those obtained 

by functionally differentiating W^i^ = —i In Z^r^ with respect to the external currents twice. Expanding in powers 
of 1/N, one can show that all the diagrams representing the corrections, as given by Eq. I|C19() . to the connected 
part of the generating functional, PFctp = —ilnZcTP, are of order or higher |53|- Therefore, one can conclude 

that the leading order contribution to Wctp is entirely given by W^ctp^ which is of order 1/A^ and from which the 
leading order contribution to all the quantum correlation functions with an even number of points can be obtained. 
Two particular examples showing how the corrections due to S'int [Sob, Acd] contribute to the two-point quantum 
correlation functions (corresponding to terms in Wctp which are quadratic in the external currents) are provided in 
Fig. n The first diagram involves vertices with three legs associated with cubic terms in the gravitational action. 
The second diagram involves a nonlocal vertex with four legs associated with quartic terms in the influence action. 
The nonlocal vertex has been represented by a loop of the matter fields because, if Feynman diagrams are introduced 
when evaluating the influence action, the terms quartic in the metric perturbations giving rise to the second diagram 
in Fig. ^ correspond to a loop of matter fields with four insertions linear in the metric perturbations (there are three 
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0(l/iV2) 



FIG. 1: Two diagrams illustrating the fact that including either the vertices for the metric perturbations (as in the first 
diagram) or terms from the influence functional evaluated beyond the Gaussian approximation (as in the second diagram) leads 
to contributions of higher order in 1/A''. In particular, the two diagrams shown here give contributions of order to the 

two-point quantum correlation function for the metric perturbations. The plain lines represent the CTP propagators for the 
matter fields on the background spacetime and the wavy lines correspond to the GTP propagators for the metric perturbations 
obtained by functionally differentiating twice with respect to the external currents. 

other contributions to the influence action involving terms quartic in the metric perturbations: one corresponds to 
a loop of matter fields with two insertions linear in the metric perturbation and a third insertion quadratic in the 
metric perturbation, a second contribution that corresponds to a loop of matter fields with two insertions quadratic 
in the metric perturbation, and a third one corresponding to a loop of matter fields with an insertion linear in the 
metric perturbation and a second insertion cubic in the metric perturbation). 

APPENDIX D: SINGULAR COINCIDENCE LIMIT FOR THE NOISE KERNEL 



The noise kernel defined by Eq. (|12|l has in general a singular coincidence limit a; — > (in fact, it is still singular 
even for x ^ x' when the two points are connected by a null geodesic), which translates into an ultraviolet divergence 
when integrating over momenta in Fourier space, as can be seen from Eq. ljB7|l . The result is, nevertheless, finite when 
X ^ x' (and they are not connected by a null geodesic). In fact, even though the noise kernel is not well defined as a 
tensor- valued function [o^ , it is well defined as a tensor- valued distribution and yields finite results when integrated 
with suitable test functions. 

Let us consider a specific example to illustrate the points addressed in this Appendix: a massless conformally 
coupled scalar field in Minkowski spacetime. The expression for the noise kernel in spacetime coordinates, which 
results from Fourier transforming Eq. ljB7|) and is a well-defined distribution, is the following 54] : 



M^,upa{x~x') CX {d^d,d'pd',)Vf 



1 



{x - x'y 



1 
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{d^d,d'd',)a,a'jn{x-x' 



/^2 



(Dl) 



where 7^/ stands for the Hadamard finite part prescription, whose precise definition can be found in Refs. |73l l74j . 
However, the contribution from the induced fiuctuations to the symmetrized two-point correlation function, which is 
given by 



— 2 

^ (Crcf • (Grct)^) {X1,X2), 



(D2) 



is not necessarily well-defined if the time integral in • involves a finite initial time U. That is because in that case 
the noise kernel is actually convoluted with Grct{xi, x'i)9{t'i — ti), which is not a good test function since it is not 
differentiable at t'^ = ti. The fact that the result for expression l|D2(l is singular for a finite initial time can actually 
be seen by using the last equality in Eq. IjDip . substituting into expression (|D2|I and integrating by parts. The 
contributions from the boundary terms at the finite times ii and ^2 are finite, at least when ti ^ t2 (the fact that 
ti t2 may be required to get a finite result is simply indicating that expression (jD2l) is also a distribution). On the 
other hand, the boundary terms that correspond to taking t'^ and t'2 both equal to ti are divergent. The fact that all 
the singular contributions can be concentrated at the initial time seems to suggest that the origin of the problem may 
be related to the initial state that was chosen. 
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We proceed now to argue that the origin of the singularities described in the previous paragraph can indeed be 
traced back to the initial state that was considered, with the metric perturbations and the matter fields completely 
uncorrelated. In order to do that, it will be useful to discuss an analogous situation for QBM models such as that 
described in Appendix^ In particular, let us consider an ohmic distribution for the environment frequencies with 
an ultraviolet cut-off A, which can be characterized by a spectral density function such as /((jj) = ^^(A — uj) or 
Lu exp{—LU / A) (the details about the particular way in which the cut-off is implemented are not important here). We 
have an expression analogous to ljD2|) for the induced fluctuations (see Eq. (jA2p and Ref. H^). If we consider the 
ground state (thermal state at zero temperature) as the initial state for the environment, the noise kernel is given 
by N{t,t') — J du I {uj)/uj cos uj{t — t'). When taking the limit A — > oo the noise kernel becomes proportional to 
Vf(l/{t'i — ^2)^)- To obtain the correlation function we integrate by parts, as described above, in the expression for 
the induced fluctuations, Eq. (|A2|I . before taking the limit A ^ 00. We obtain again a boundary term at the initial 
time which diverges as we finally let A go to infinity, and we end up with an infinite result for the correlation function. 
On the other hand, one can show that the result for the correlation function of the ground state of the whole system 
(system plus environment), including the system-environment interaction, is finite (the use of Euclidean path integrals 
is particularly convenient in this respect) |75j |. This constitutes a clear example of the fact that initial states in which 
the system and the environment are suitably correlated give rise to well-defined correlation functions. 

Alternatively, when taking a completely uncorrelated initial state, one can still get a finite result for the correlation 
function by smoothly switching on the system-environment interaction so that the boundary term at the initial time 
which results from the integration by parts and becomes divergent in the limit A — > 00 actually vanishes. This reveals 
again that the origin of the singularity for the correlation function arises because the highest frequency modes of 
the environment become correlated with the system in a timescale of the order of A^^. Such a fact is supported by 
the existence of a jolt with a characteristic timescale A"-'^ in the diffusion coefficients of the master equation which 
becomes singular when A 00, as was found in Ref. [t^. In fact, one can show that those states in which the high 
frequency modes of the environment and the system are uncorrelated are unphysical when the environment contains 
an infinite number of modes with arbitrarily high frequencies since their energy becomes infinite as A — s- cx). 

Return now to the gravitational case. There are some situations, such as the effect of stress tensor fluctuations 
on the propagation of null geodesies, in which the appropriate way to deal with the singular coincidence limit of the 
noise kernel is by integrating over some smearing function Q (in general smearing just along the spatial directions is 
not enough: smearing in time is needed to get a finite result). On the other hand, when computing the correlation 
functions for the metric perturbations, the noise kernel naturally appears integrated with the retarded propagator. As 
explained above, the problem still persists at the initial time, which reflects the unphysical character of the completely 
uncorrelated initial state that was employed. Similar to QBM models, a well-defined result for the correlation functions 
can be obtained by considering a properly correlated initial state, such as that resulting from the use of Euclidean 
path integrals that are then analytically continued to Lorentzian time |6^ . Roughly speaking, this would imply 
the existence of an additional term in Eqs. (|14ll and IjClSI) due to the existence of correlations between the initial 
conditions for the solutions of the Langevin equation and the stochastic source, which reflect the initial correlations 
between the system and the environment. 

Alternatively, one can still make sense of the results obtained from assuming an uncorrelated initial state by smoothly 
switching on the interaction between the metric perturbations and the matter fields so that the high frequency modes 
can get correlated with the system. However, in contrast to the QBM case, we have to be careful with switching on 
the interaction during a finite period of time since that would imply that the source of the Einstein-Langevin equation 
is not conserved and would be in conflict with the Bianchi identity, which guarantees the integrability of the equation. 
Therefore, the interaction should be turned on adiabatically and asymptotically past initial conditions should be 
considered. In fact, in Sec. IIVI where we assumed asymptotic initial conditions and worked mostly in Fourier space, 
a finite result for the correlation function was obtained without the need for explicitly switching on the interaction 
adiabatically. There are, however, situations (for instance, in cosmology) in which asymptotic initial conditions are 
not adequate. An alternative procedure should be considered in those cases. 

APPENDIX E: RUNAWAY SOLUTIONS AND METHODS TO DEAL WITH THEM 

In this Appendix we will briefly discuss the existence of runaway solutions in SCG (solutions which grow without 
bound in timescales comparable to the Planck time), their counterparts at the quantum level, and how their connection 
can be understood in the context of stochastic gravity. We will also discuss the existing prescriptions for dealing with 
this kind of unstable solutions. 
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1. Runaway solutions in semiclassical gravity 



Let us start by considering the linearized semiclassical Einstein equation around the Minkowski spacetime. The 
solutions for the case of a massless scalar field were first discussed in Ref. |0 and an exhaustive description can 
be found in Appendix A of Ref. [^. Taking Eq. (Tj) and using a decomposition for the linearized Einstein tensor 
analogous to that introduced in Sec. IIV CI for the metric perturbation, the vectorial part is found to vanish [o^. 
whereas the scalar and tensorial contributions satisfy the following equations: 

(Fi (p) + 3p^F2 ip)) Gi'} ip) = 0, (El) 
Fi{p)Gi'}^^\p)^0. (E2) 

where Fi{p) and F2{p) are given by Eqs. Q-CDl), and G/jj and G/ij denote, respectively, the scalar and tensorial 
parts of the linearized Einstein tensor. In order to illustrate how the runaway solutions arise, we will consider the 
particular example of a massless and conformally coupled scalar field (see Ref. |22| for the massless case with arbitrary 
coupling and Refs. 0,01 for the general massive case). The previous equations become then 

(l + 12«:^p2)(5W(S)(^)^0^ (23) 
lim^ (l + (9607r2)-i«;p2 ( ) ) '■^\p) = 0. (E4) 

In addition to the obvious solution g\}u ^^"^ (p) = (the only solution when /3 = 0) , when /? > the solutions for 
the scalar component exhibit an oscillatory behavior in spacetime coordinates which corresponds to a massive scalar 
field with = (12k|^|)~^; for ^ < the solutions correspond to a tachyonic field with — ~{12k\(3\)~^: in 
spacetime coordinates they exhibit an exponential behavior in time -growing or decreasing- for wavelengths larger 
than 47r(3K|/3|)^/^, and an oscillatory behavior for wavelengths smaller than 47r(3K|^|)^/^. On the other hand, the 
solution G^i}}''^\p) = is completely trivial since any scalar metric perturbation h^y{p) giving rise to a vanishing 
linearized Einstein tensor can be eliminated by a gauge transformation as explained in Sec. IIV CI 

As for the tensorial component, when < /^crit = /j7^(1207r)^/^e^ (or A > Adit = in the notation of Ref. ;2lJ) 
the first factor in Eq. (jE4|) vanishes for four complex values of p^ of the form and ±w* , where w is some complex 
value, as illustrated in Fig. [3 We will consider here the case in which /i < //crit; a detailed description of the situation 
for p > /icrit can be found in Appendix A of Ref. [2ll | . The two zeros on the upper half of the complex plane correspond 
to solutions in spacetime coordinates exponentially growing in time, whereas the two on the lower half correspond 
to solutions exponentially decreasing in time. Strictly speaking, these solutions only exist in spacetime coordinates, 
since their Fourier transform is not well defined. They are commonly referred to as runaway solutions and for /i ~ 
they grow exponentially in timescales comparable to the Planck time. 



2. Quantum mechanical systems with higher order time derivatives 

Before proceeding to discuss the situation in stochastic gravity, it is interesting to make a few remarks about the 
quantization of higher derivative theories and the counterparts of the previous classical instabilities in the quantum 
context. Let us consider first a free theory with a structure analogous to that of linearized semiclassical gravity around 
Minkowski spacetime without including the nonlocal terms. It is characterized by the following Lagrangian, which 
corresponds to a harmonic oscillator with a higher derivative term: 

L{q,q,q) = ^f + ^e-^^'q^ (E5) 

To begin with, one can consider a generalization of the usual canonical formalism introduced by Ostrogradski to deal 
with theories involving higher order derivatives (see, for instance Refs. 19, 77]). The theory can then be quantized 
following the standard canonical quantization rules. The corresponding Wigner function (or Wigner functional if 
a field theory were considered) can also be introduced. The pathological character of the theory becomes clear by 
diagonalizing the Hamiltonian and realizing that the result corresponds to two independent harmonic oscillators, but 
with one of them having a negative sign in the kinetic term. For r < the potential term of the harmonic oscillator 
with the negative kinetic term is also negative and the classical solutions do not exhibit instabilities. However, in 
any case the configurations for the harmonic oscillator with the negative kinetic term can have negative energies 
arbitrarily large in absolute value. Moreover, the frequency for that oscillator is proportional to r~^/^ and diverges 
as T — > 0. At the quantum level, such a theory also gives rise to negative eigenvalues of the Hamiltonian arbitrarily 
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large in absolute value, but can be alternatively formulated in terms of a Hamiltonian without negative energies by 
introducing states with negative norm (commonly referred to as ghosts) . This fact is often argued in a qualitative 
way by pointing out that the propagator of the theory in Fourier space is proportional to 

^TT—T- (E6) 



It should also be mentioned that Hawking and Hertog have suggested a prescription for dealing with that kind of 
theories which is based on imposing well-defined boundary conditions in Euclidean time and then Wick rotating back 
to Lorentzian time. The results have then a nonsingular limit t ^ 0, so that when the higher order derivative term 
in the Lagrangian is small, one essentially recovers the results of the second order theory 65]. 

Even though there is a range of parameters (t < 0) in which the free theory described above does not exhibit 
instabilities, they arise when a nonlinear self-interaction term is added to the Lagrangian. The reason is that the two 
Hamiltonian contributions corresponding to a couple of harmonic oscillators, one with a negative energy spectrum 
and the other with a positive one, can have a stable evolution as long as they are decoupled. However, adding an 
interaction term couples them in such a way that one can acquire negative energies arbitrarily large in absolute value 
while the other gains large positive energies, which is the source of instability. In general this is reflected in the 
structure of the propagator as a shift of the poles on the real axis to the complex plane. Hawking et al. have argued 
that well-behaved results can still be obtained by imposing boundary conditions which discard solutions which grow 
unboundedly in time (64,, ,,65j . Those conditions can be implemented by a suitable choice of the integration contour on 
the complex plane when computing the inverse Fourier transform of the propagator, but causality is violated at small 
timescales (we will come back to this point below). Another possibility, when the parameter r is small, is to make use 
of an order reduction procedure ^2L .22,. 78]. which consists of differentiating the equation of motion with respect to 
time, substituting back into the original equation and discarding the terms of higher order in r. This procedure can 
be iterated as many times as necessary to get a second order equation valid up to the corresponding order of r. The 
usual canonical formalism associated with the second order equation of motion can then be employed to evolve the 
Wigner function. It should be stressed that, although we have considered a simple model as an illustrating example, 
the previous methods have been applied to more involved situations, including SCG plll22| and quantum cosmology 

il. 



3. Runaway solutions in stochastic gravity 



Let us now address the case of stochastic gravity and see how the instabilities in SCG and the difhculties in quantizing 
theories with higher order derivatives are related. First of all, we recall that in Appendix [CI the counterterms quadratic 
in the curvature were ignored and it was implicitly assumed that the Einstein-Langevin equation was a second order 
intcgro-differential equation whose initial conditions were completely determined by specifying the metric perturbation 
and its normal derivative on the initial Cauchy hypersurface. If the counterterms quadratic in the curvature, which 
give rise to higher order derivative terms, are also taken into account, the generalized canonical formalism referred 
to above and the corresponding Wigner functional should be used. In fact, due to the singular behavior of the 
nonlocal part of the dissipation kernel at the initial time, specifying initial conditions at a finite initial time is an 
even more delicate matter. In any case, since we have to consider asymptotic initial conditions to deal with the 
singular coincidence limit of the noise kernel, as explained in Appendix^ we do not need to be concerned about the 
problems associated with finite initial times. Runaway solutions, however, still exist and some method to deal with 
them is required. In particular, when computing two-point correlation functions in the context of stochastic gravity, 
the existence of runaway solutions has implications for both the intrinsic and the induced contributions. 

One possible method for dealing with the existence of runaway solutions is the order reduction prescription. As 
explained above, the method is based on treating perturbatively the terms involving higher order derivatives, differ- 
entiating the equation under consideration and substituting back the higher derivative terms in the original equation 
keeping only terms up to the required order in the perturbative parameter. In the case of the semiclassical Einstein 
equation, the perturbative parameter employed is h or, equivalently, the square of the Planck length = k/Stt. If we 
consider the semiclassical Einstein equation for linear metric perturbations around Minkowski spacetime and differ- 
entiate twice with respect to the background covariant derivative, it becomes clear that the second order derivatives 
of the Einstein tensor are of order k. Substituting back into the original equation, we get the following equation up 
to order k^: 

G^2[9 + h] = Q + 0{K'), (E7) 

where no effects from the vacuum polarization of the quantum matter fields are left. Since the linearized semiclassical 
Einstein equation coincides with the homogeneous part of the Einstein-Langevin equation, Eq. (jE7|) governs the 
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contribution of the intrinsic fluctuations to the quantum correlation function, which coincides with that of free 
gravitons. Similarly, when making use of the order reduction prescription, the Einstein-Langevin equation becomes 

G'^^^[g + h] = KU + 0{>^), (E8) 

where the stochastic source, whose correlation function only depends on the background metric and hence does not 
involve higher order derivatives of the metric perturbation, is not affected by the order reduction procedure. Therefore, 
in contrast to the intrinsic fluctuations, there will still be a nontrivial contribution to the induced fluctuations due to 
the polarization of the quantum matter fields, but no contribution from the dissipation kernel is left in the Einstein- 
Langevin equation. Since all the terms involving higher order derivatives, which were associated with the dissipation 
kernel, have been discarded, an ordinary Wigner functional can be introduced without any need to consider generalized 
Ostrogradski momenta. Furthermore, the absence of the dissipation kernel also allows the possibility of specifying 
initial conditions at a finite initial time as far as the homogeneous solutions (relevant for the computation of the intrinsic 
fluctuations) and the retarded propagator are concerned. Nevertheless, one is still forced to consider asymptotic initial 
conditions in order to get a finite result for the induced fluctuations due to the singular coincidence limit of the noise 
kernel, as explained in the previous Appendix. 

Hawking et al. have proposed an alternative procedure for dealing with the runaway solutions 0, Ig^. Their 
method is based on imposing final boundary conditions which discard those solutions that grow unboundedly in time. 
Let us first see how their approach can be applied to the computation of the intrinsic fluctuations by considering 
the particular case of a massless and conformally coupled scalar field. From Eq. (|E4|I and Fig. [21 one can see that, 

in addition to the solution gJ^] = 0, the solutions of the tensorial part grow or decrease exponentially in time. The 
exponentially growing solutions are discarded when the final boundary condition is imposed, and the contributions 
from the exponentially decreasing ones also vanish if regular initial conditions are specified at an asymptotic initial 
time. On the other hand, from Eq. I)E3|I one can see that the situation is analogous for the solutions of the scalar part 
when f3 < 0. For f3 > the solutions are oscillatory and, hence, are not discarded when the final boundary condition 
is imposed (in contrast to the situation where the order reduction prescription is used). 

Let us now apply the previous approach to the computation of the induced fluctuations. When considering asymp- 
totic initial conditions, the relevant propagator for expressing the linearized Einstein tensor in terms of the stochastic 
source, can be obtained by inverting F^^apiv) in Eq. |(7J). The resulting propagator, D^^a[}{p), exhibits a number of 
poles in the complex plane, as illustrated in Fig. |21 The expression for the retarded propagator in spacetime coordi- 
nates corresponds to choosing the integration path represented by the dashed line in Fig. [3 when Fourier transforming 
back from momentum space. It exhibits the appropriate causal behavior: D^^^apix — y) ~ for tx < ty, as can be 
seen by closing the integration contour on the upper half of the complex plane. However, for > ty it increases 
exponentially in time due to the contributions from the two poles on the upper half of the complex plane when closing 
the path on the lower half. Imposing the final boundary conditions which discard solutions growing unboundedly 
in time is equivalent to taking a different integration path: that represented by a solid line in Fig. |21 The resulting 
propagator does not exhibit exponential instabilities, but gives rise to causality violations since Df^^ai3{x — y) ^ for 
tx < ty (the characteristic timescale of these causality violations is of order \/Nlp). This propagator is the only one 
which has a well-defined Fourier transform. It was employed in Ref . , where it was argued that any other propaga- 
tor should yield an equivalent result for the correlation function obtained by solving the Einstein-Langevin equation. 
This argument is certainly true for propagators with a well-defined Fourier transform. However, the existence of poles 
off the real axis gives rise to propagators in spacetime coordinates (they do not have a well-defined Fourier transform 
because of the exponentially growing or decreasing contributions) which yield inequivalent results for the correlation 
function. Since this choice for the propagator was made, the results obtained in Ref. [sJl correspond to those that 
would follow when employing the procedure proposed by Hawking et al. In fact. Hawking et al. applied their method 
to quantum propagators, but, as we have described, it can also be used when solving the semiclassical Einstein equa- 
tion and the Einstein-Langevin equation. The stochastic correlation functions obtained are then equivalent to the 
quantum correlation functions (CTP propagators) which would result from the application of the prescription. 

4. Estimates of radiative corrections for a single matter field and a large number of them 

SCG is expected to provide reliable results as long as the characteristic lengthscales under consideration are much 
larger than the Planck length Ip . This can be qualitatively argued by estimating the magnitude of the different 
contributions to the effective action (considering the relevant Feynman diagrams and using dimensional arguments): 
the Einstein-Hilbert term and the radiative quantum corrections. The Einstein-Hilbert term is of order Ip'^R (the 
characteristic curvature R is simply given by L~^, where L is the characteristic lengthscale of our problem); the 
vacuum polarization terms involving loops of matter fields are of order R^; and higher loop corrections involving 




FIG. 2: Representation in the complex plane of the values of p for which the coefflcient of the semiclassical Einstein equation 
for the tensorial components of the Einstein tensor in Fourier space vanishes. The case fj, < ^crit and a particular value of 
were considered, but the qualitative structure will remain the same for any other value of p^. The plot also corresponds to the 
structure of the poles of the tensorial part of the propagator D^^ci/3(p). The solid line corresponds to the integration contour 
when Fourier transfor ming back to spacetime coordinates which follows from the prescription proposed by Hawking et al. 
and was chosen in Ref. ||54| . Changing this contour as indicated by the dashed lines, one obtains a strictly retarded propagator, 
but it exhibits exponential instabilities for large positive time differences associated with the two poles on the upper half of the 
complex plane. 



internal graviton propagators are of order IpR^ or higher. Thus, we see that the higher order corrections not included 
in SCG are neghgible provided that L ^ Ip. In that regime, however, the vacuum polarization terms only yield a small 
correction to the Einstein-Hilbert term and any classical gravitational source which were present. The justification 
of the order reduction prescription is actually based on this fact. Therefore, significant effects from the vacuum 
polarization of the matter fields are only expected when their small corrections accumulate in time, as would be the 
case, for instance, for an evaporating macroscopic black hole all the way before reaching Planckian scales. 

The previous estimates for the different terms in the effective action change in a remarkable way when a large 
number of fields, N, is considered |95j| . The vacuum polarization terms involving loops of matter become of order 
NR^ and, similarly, the higher loop corrections involving internal graviton propagators are of order Nl^R^ or higher 
(the contributions corresponding to one and two graviton loops are, respectively, of order and IpR^, but are 
negligible as compared to those from matter loops when N is large). There is then a regime in which the vacuum 
polarization of the matter fields and the Einstein-Hilbert term are comparable when L ^ ^/Nlp. On the other hand, 
the higher loop corrections will still be much smaller if L ^ Ip. Both conditions are compatible provided that the 
number of fields, N , is very large. This is, in fact, the kind of situation considered in trace anomaly driven inflationary 
models '64] , such as that originally proposed by Starobinsky '6^ , where the exponential inflation is driven by a large 
number of massless conformal fields. The order reduction prescription would completely discard the effect from the 
vacuum polarization of the matter fields even though it is comparable to the Einstein-Hilbert term. In contrast, the 
procedure proposed by Hawking et al. keeps the contribution from the matter fields. 

We conclude this Appendix by mentioning that it has been pointed out that a similar kind of instability, which 
is closel y co nnected to the existence of the Landau pole, is also present in scalar QED (as well as ordinary QED) 
[TtL llsL l3^. Nevertheless, a number of nonperturbative studies on the evolution of the expectation value of the 
electromagnetic fleld using a large N expansion have been carried out. In fact, it was suggested in Ref. ji^ that by 
introducing a flnite 3-momentum cut-off and considering a running coupling constant small enough at low energies, 
the problem with the Landau pole could be circumvented (at least from a practical point of view). Yet it seems 
unlikely that a similar procedure could work for the gravitational case due to the existence of higher derivatives. 
Moreover, introducing a 3-momentum cut-off would break general covariance and that would pose serious difficulties 
when implementing a consistent and natural renormalization scheme in general curved spacetimes. 
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[88] In those References dealing with stochastic gravity this term is usually called the dissipation term by analogy with QBM 
models. 

[89] The existence of such a decay still holds for the overdamped case provided that fJrcn 7^ 0, otherwise there is a constant 

contribution that does not decay in time. 
[90] Strictly speaking, the term dissipation kernel commonly refers to the antisymmetric part of the kernel —Hicn- Making 

an abuse of language, we will employ this term to refer to the whole kernel — 2i/ren plus the local and symmetric kernel 

-2Mron. 

[91] In general, one should be careful with the contributions from the timelike boundaries as well as the edges connecting the 
spacelike and timelike boundaries 79). Here we will assume that the timelike boundaries are infinitely far away and the 
value of the metric perturbations decays at large distances so that only the contributions from the spacelike boundaries 
are relevant. 

[92] Throughout this Section we will neglect the contribution to the momentum from the counterterms; see Appendix IeI for 
further discussion on this point. 
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[93] The noise kernel is in general a bitensor, but due to the homogeneity of Minkowski spacetime and the triviaUty of the 
connection (and the corresponding parallel transport), the noise kernel becomes in that case a simple tensorial field which 
depends on {x — x')'^ . 

[94] More precisely, decomposing the metric perturbation into scalar, vectorial and tensorial parts, as done in Sec. IIV CI and 
computing the linearized Einstein tensor, one gets a vanishing result for the vectorial part of the metric perturbation; the 
scalar and tensorial components of the metric perturbation give rise, respectively, to the scalar and tensorial components 
of the linearized Einstein tensor. 

[95] The actual physical Planck length Ip is considered, not the rescaled one, ^ k/Stt, which is related to Ip by Snip = k — R/N . 



